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ABSTRACT 


Arakawa  and  Schubert  (1974)  developed  a  cumulus  parameterization  scheme  in 
a  framework  chat  conceptually  divides  Che  mutual  interaction  between  Che 
cumulus  convection  and  large-scale  disturbance  into  the  categories  of  large- 
scale  budget  requirements  and  Che  quasi-equilibrium  assumption  of  cloud  work 
function.  We  have  applied  the  Arakawa-Schubert  scheme  through  a  semi¬ 
prognostic  approach  to  two  different  data  sets:  one  is  for  an  intense  tropical 
cloud  band  event;  the  other  is  for  tropical  composite  easterly  wave 
disturbances.  Both  were  observed  in  GARP  (Global  Atmospheric  Research  Program) 
Atlantic  Tropical  Experiment  (GATE).  The  cloud  heating  and  drying  effects 
predicted  by  the  Arakawa-Schubert  scheme  are  found  to  agree  rather  well  with 
Che  observations.  However,  it  is  also  found  that  Che  Arakawa-Schubert  scheme 
underestimates  both  condensation  and  evaporation  rates  substantially  when 
compared  with  the  cumulus  ensemble  model  results  (Soong  and  Tao,  1980;  Tao, 
1983).  An  inclusion  of  the  downdraft  effects,  as  formulated  by  Johnson  (1976), 
appears  to  alleviate  this  deficiency. 

In  order  to  examine  how  the  Arakawa-Schubert  scheme  works  in  a  fully 
prognostic  problem,  a  simulation  of  the  evolution  and  structure  of  the  tropical 
cloud  band  mentioned  above  under  the  influence  of  an  imposed  large-scale,  low- 
level  lifting  was  made  by  using  a  two-dimensional  hydrostatic  model  with  the 
inclusion  of  the  Arakawa-Schubert  scheme.  Basically,  the  model  result 
indicates  that  the  mesoscale  convective  system  is  driven  by  the  excess  of  the 
convective  heating  derived  from  the  Arakawa-Schubert  scheme  over  the  adiabatic 
cooling  due  to  the  imposed  large-scale  lifting  and  induced  mesoscale  upward 
motion.  However,  as  the  convective  system  develops,  the  adiabatic  warming  due 
to  the  subsidence  outside  the  cloud  system  gradually  accumulates  into  a 
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secondary  temperature  anomaly  which  subsequently  reduces  the  original 
temperature  contrast  and  inhibits  the  further  development  of  the  convective 
system.  A  24-hour  integration  shows  that  the  model  is  capable  of  simulating 
many  Important  features,  such  as  the  life  cycle,  the  intensity  of  circulation, 
and  rainfall  rates. 


CHAPTER  I 
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t. 
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INTRODUCTION 


The  importance  of  the  vertical  transports  of  heat,  moisture,  vorticlty, 
and  momentum  by  cumulus  clouds  has  long  been  recognized  not  only  in  the 
numerical  predictions  of  weather  and  climate  but  also  In  the  simulation  studies 
on  mesoscale  circulations  such  as  tropical  convective  systems,  mid-latitude 
f rontogenesis ,  and  the  sea-land  breeze.  However,  it  is  impractical  to  resolve 
those  Individual  cumulus  clouds  because  of  the  distinct  scale  difference 
between  them  and  larger-scale  circulations.  In  the  past  two  decades,  there 
have  been  numerous  attempts  to  relate  subgrid-scale  cumulus  convection  to 
large-scale  fields  by  specifying  the  cumulus  activity  as  a  function  of  large- 
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scale  variables,  which  is  now  commonly  referred  to  as  cumulus 
parameterization".  The  goal  of  cumulus  parameterization  is  to  predict  the 
change  of  large-scale  disturbances  due  to  cumulus  convection  by  describing  not 
each  of  the  many  individual  clouds,  but  only  their  collective  effects. 

The  first  attempts  to  parameterize  cumulus -convection  have  been  known  as 
Conditional  Instability  of  the  Second  Kind  (CISK)  by  Charney  and  Eliassen 
(1964)  and  Ooyama  (1964).  The  concept  of  CISK,  which  was  particularly  designed 
to  explain  the  growth  of  tropical  cyclones,  can  be  summarized  as  follows:  the 
release  of  latent  heat  generated  by  cumulus  cloud  drives  the  cyclone  and,  in 
turn,  the  cyclone  provides  the  moisture  which  maintains  the  cumulus  clouds. 
Later,  many  studies  extended  the  concept  of  CISK  in  investigating  tropical  wave 


disturbances  through  normal-mode  linear  analysis,  which  is  commonly  called  the 
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Wave-CISK  theory.  The  study  by  Lindzen  (1974)  is  representative  of  the  state 
of  the  art.  However,  Ooyama  (1982)  recently  made  a  thorough  review  of  Che 
validity  and  limitations  of  CISK  and  pointed  out  that  many  arguments  about 
"preferred  scales"  in  the  linear  theory  were  irrelevant  to  the  tropical 
cyclones,  especially  to  their  mature  stages  where  the  nonlinear  processes  are 
dominant . 

Since  the  emergence  of  CISK  theory,  a  great  number  of  methods  or  "schemes" 
dealing  with  the  interaction  of  cumulus  convection  with  the  large-scale 
disturbances  have  been  developed.  Each  involves  different  specific 
formulations  based  on  the  model's  needs  and  the  author's  intuitive 
idealizations  of  cumulus  convection.  To  date,  there  are  three  major  kinds  of 
cumulus  parameterization  schemes  which  have  been  more  generally  employed  in 
large-scale  and  mesoscaie  numerical  models.  These  are:  (1)  Convective 
adjustment  scheme  (Manabe  £]t  aj^. ,  1965);  (2)  Kuo's  Scheme  (Kuo,  1965  ,  1974); 
and  (3)  the  Arakawa-Schubert  scheme  (Arakawa  and  Schubert,  1974). 

The  concept  of  the  convective  adjustment  scheme  is  rather  simple  and 
implicit.  Whenever  the  atmosphere  is  super-saturated  and  the  temperature 
stratification  is  unstable  with  respect  to  the  moist  adiabatic  lapse  rate,  the 
scheme  assumes  that  cumulus  clouds  form  and  redistribute  the  temperature  and 
moisture  profiles  so  that  the  total  moist  enthalpy  is  conserved  and  the  lapse 
rate  is  moist  adiabatic.  This  scheme  has  been  extensively  used  in  the  general 
circulation  model  (CCM)  at  the  Geophysical  Fluid  Dynamics  Laboratory  (GFDL) . 
This  method  was  also  adopted  by  the  National  Center  for  Atmospheric  Research 
(NCAR)  for  its  Community  Climate  Model  (CCM).  In  addition,  Kurihara  (1973)  and 


Kurihara  and  Tuleya  (1974)  applied  a  modified  version  of  this  adjustment  scheme 
to  hurricane  modeling  studies. 

Kuo's  scheme  (1965)  parameterizes  the  release  of  latent  heat  by  relating 
the  convective  activity  to  the  large-scale  moisture  convergence  into  a 
conditionally  unstable  atmospheric  column.  The  adjustment  of  the  large-scale 
temperature  and  moisture  fields  due  to  clouds  can  be  determined  from  the  large- 
scale  horizontal  convergence  of  moisture  and  the  profiles  of  excess  temperature 
and  mixing  ratio  inside  the  cloud.  Kuo  (1974)  modified  his  scheme  by 
separating  the  total  convergence  of  moisture  into  two  parts:  one  (denoted  by 
b  X  Mj.)  for  increasing  the  humidity  of  the  air;  the  other  (denoted  by 
(1-b)  X  M^)  for  warming  the  air  by  Che  release  of  latent  heat.  There  have  been 
quite  a  few  studies  which  attempted  to  estimate  the  value  of  b  (Cho,  1976; 
Krishnamurti  et_  aj^. ,  1980,  etc.),  but  no  satisfactory  theoretical  determination 
for  b  is  available  yet.  Kuo's  1965  scheme  has  been  extensively  used  at  the 
National  Meteorological  Center  (NMC)  for  various  operational  and  research 
models.  The  European  Center  for  Medium  Range  Forecasts  is  currently  using 
Kuo's  1974  Scheme.  Theoretically  speaking,  both  the  convective  adjustment 
scheme  and  Kuo's  scheme  are  highly  empirical,  neglecting  many  of  the  physical 
processes  involved  in  the  mutual  interaction  between  clouds  and  Che 
environment.  Consequently,  these  two  schemes  are  not  able  to  provide  detailed 
information  about  the  properties  of  cumulus  ensembles  such  as  cloud  mass  flux, 
liquid  water  detrainment,  condensation  and  evaporation  rates  and  vertical  eddy 
transports  of  heat  and  moisture.  Hence,  from  the  academic  point  of  view,  the 
convective  adjustment  scheme  and  Kuo's  scheme  do  not  serve  the  whole  purpose  of 
cumulus  parameterization. 


Ooyama  (1971)  first  proposed  a  theory  taking  into  account  the  coexistence 
of  a  spectrum  of  clouds  which  were  represented  by  independent  entraining 
buoyant  elements  dispatched  from  the  mixed  layer.  However,  this  theory  was  not 
closed  since  the  determination  of  the  dispatcher  function  was  left  to  future 
consideration.  Based  on  the  concept  of  Ooyama's  spectral  clouds  and  the 
earlier  works  of  Arakawa  (1969,  1972),  a  theory  of  the  interaction  between  a 
cumulus  ensemble  and  the  large-scale  environment  has  been  developed  by  Arakawa 
and  Schubert  (1974).  In  essence,  the  Arakawa- Schubert  scheme  divides  the 
cumulus  ensemble  into  one-dimensional,  steady-state,  entraining  sub-ensemble 
clouds;  and,  more  importantly,  the  scheme  is  closed  by  a  physically  well- 
defined  hypothesis.  This  scheme  can  be  partitioned  into  two  parts:  (1)  the 
budget  requirements  for  the  large-scale  disturbances,  which  describe  that  the 
cumulus  convection  affects  the  temperature  and  moisture  fields  through  cumulus- 
induced  subsidence,  detrainment  of  saturated  air,  and  evaporation  of  detrained 
liquid  water;  (2)  the  assumption  of  quasi-equilibrium  for  the  cloud  work 
function,  which  describes  how  the  clouds  stabilize  the  atmosphere  while  the 
large-scale  motion  generates  moist  convective  instability. 

The  "quasi-equilibrium  assumption"  permits  the  determination  of  the  cloud 
mass  flux  at  cloud  base  (a  unique  unknown  variable  in  tlie  scheme)  so  that  the 
large-scale  budget  equations  can  be  used  in  a  prognostic  sense.  Meanwhile,  the 
vertical  distributions  of  cloud  mass  flux,  cloud  detrainment,  condensation  and 
evaporation  rates,  and  eddy  transports  of  heat  and  moisture  can  be  also 
determined.  Besides,  the  quasi-equilibrium  assumption  involves  the 
significance  of  the  interaction  between  clouds.  Thus,  as  far  as  the  physical 
processes  are  concerned,  the  Arakawa-Schubert  scheme  has  been  considered  Che 
most  complete  cumulus  parameterization  scheme.  This  report  will  focus  on  the 


Arakawa-Schubert  scheme  (hereafter  referred  to  as  the  A-S  scheme)  not  only 
because  it  contains  more  sophisticated  physical  processes  but  also  because  it 


has  been  explored  less  frequently  in  the  past  when  compared  with  the  convective 
adjustment  scheme  and  Kuo's  scheme. 

Investigations  with  the  use  of  the  A-S  scheme  in  linear  Wave-CISK  models 
have  been  conducted  by  Stuart  (1976),  Shukla  (1978),  and  Crura  and  Stevens 
(1983).  Since  the  interaction  between  clouds  and  large-scale  disturbances  is 
inherently  a  nonlinear  process  and  since  the  linear  normal-mode  concept  can 
cause  negative  heating  which  is  physically  unrealistic,  we  believe  that 
additional  research  on  the  role  of  the  A-S  scheme  in  the  interaction  between 
clouds  and  the  environment  should  be  carried  out  with  a  numerical  prognostic 
model.  Miyakoda  and  Sirutis  (1977)  compared  the  convective  adjustment  scheme 
in  the  GFDL  GCM  with  the  A-S  scheme  in  the  UCLA  (.CM.  Basically,  they  found 
that  the  A-S  scheme  can  produce  deeper  penetrative  convective  clouds.  Also, 
the  distribution  of  rainfall  relative  to  the  trough  and  ridge  of  tropical 
easterly  waves  is  different  between  the  two  schemes.  However,  in  this 
experiment,  not  only  the  cumulus  parameterization  schemes  but  also 
parameterization  schemes  for  other  processes  are  different.  Consequently,  it 
is  not  entirely  possible  to  isolate  the  differences  due  to  the  different 
cumulus  parameterizations .  Krishnamurti  a^.  (1980)  evaluated  the  rainfall 
rates  produced  by  several  different  schemes  with  the  CATE  Phase  III  data  in  a 
semi-prognostic  approach  (i.e.,  one  time-step  integration).  They  concluded 
that  Kuo's  1974  scheme  and  the  A-S  scheme  provided  the  least  root-mean-square 
error  between  the  calculated  and  observed  estimates.  However,  neither  of  these 
comparative  studies  enlightened  readers  on  the  intrinsic  nature  of  the  schemes 


themselves.  In  particular,  the  physical  reasons  behind  the  model  results  or 


the  skill  scores  derived  from  each  scheme  were  not  fully  examined.  Recently,  a 
series  of  extensive  studies  evaluating  the  quasi-equilibrium  assumption  and 
testing  the  A-S  scheme  was  conducted  by  Lord  (1978,  1982),  Lord  and  Arakawa 
(1980),  and  Lord  £t  (1982)  also  on  the  basis  of  a  semi-prognostic  approach 
using  the  GATE  (GARP  (Global  Atmospheric  Research  Program)  Atlantic  Tropical 
Experiment)  Phase  III  data.  The  predicted  rainfall  rates  by  Lord  are  generally 
in  good  agreement  with  observations. 


By  "semi-prognostic"  we  mean  that  an  instantaneous  evaluation  of  tendency 
terras  produced  by  the  cumulus  parameterization  at  a  single  grid  point  in 
response  to  a  large-scale  forcing  is  made  and  that  the  computation  is 
terminated  after  only  one  time-step.  The  advantage  of  using  the 
semi-prognostic  approach  to  test  the  cumulus  parameterization  scheme  is  that  we 
can  make  a  prediction  without  the  use  of  a  large-scale  prediction  model. 
Moreover,  since  this  is  only  a  one  time-step  integration,  the  predicted  results 
are  free  from  numerical  nonlinear  errors  and  complicated  physical  feedback 
effects  resulting  from  the  interactions  with  the  boundary  layer  and  with 
radiation  physics.  We  have  also  tested  the  A-S  scheme  semi-prognosticaily  with 
a  computat iona 1  algorithm  different  from  Lord's,  using  data  gathered  from  an 
intense  cloud  band  (cloud  cluster)  that  developed  over  the  GATE  A/B  scale  ship 
array  on  12  August  1974.  A  semi-prognostic  study  using  the  composite  GATE 
Phase  II  and  III  data  has  also  been  made.  In  both  of  these  studies,  the  cloud 
heating  and  drying  effects  predicted  by  the  A-S  scheme  were  found  to  agree  well 
with  observations.  Investigations  have  been  made  not  only  on  the  net  cloud 
heating  and  drying  effects  but  also  on  their  associated  condensation, 
evaporation,  and  vertical  transports  of  heat  and  moisture.  Furthermore,  we 


have  modified  the  A-S  scheme  by  including  the  convective-scale  downdrafts  to 


f  ^ 
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explore  how  the  cloud  properties  such  as  condensation,  evaporation,  and  the 
convergence  of  heat  and  moisture  fluxes  can  be  altered  due  to  the  downdraft 
effects.  For  comparison  purposes,  the  results  derived  from  a  cumulus  ensemble 
model  on  the  cloudband  case  by  Soong  and  Tao  (1980,  1984)  were  used. 


Even  though  the  applicability  of  the  "quasl-equlllbrlum  assumption"  to  the 
time  and  space  scales  of  tropical  convective  systems  Is  still  under  debate. 
Hack  and  Schubert  (1980)  and  Hack  et  al.  (1984)  applied  the  A-S  scheme  to  an 
axlsymmetrlc  hurricane  model  and  demonstrated  chat  the  scheme  has  Che  Inherent 
potential  to  establish  a  better  understanding  of  Che  physical  nature  of  the 
Interaction  between  organized  cumulus  convection  and  the  large-scale  fields. 
Based  on  the  encouraging  results  derived  from  our  seml-prognosClc  studies 
mentioned  above,  we  subsequently  made  another  test  on  the  A-S  scheme  using  a 
fully  prognostic  mesoscale  model.  In  order  to  examine  the  evolution  and  the 
structure  of  a  tropical  cloud  cluster  under  the  int iuence  of  a  low-level, 
large-scale  lifting,  a  two-dimensional  hydrostatic  model  was  developed  with  the 
Inclusion  of  Che  A-S  scheme.  The  model  was  applied  to  the  same  cloudband  event 
mentioned  previously.  A  time-independent ,  large-scale,  low-level  forcing  was 
imposed  on  the  model  throughout  the  whole  Integration.  The  model  results 
Indicate  chat  the  mesoscale  convective  system  is  driven  primarily  by  the 
horizontal  temperature  gradient  resuitsing  from  Che  latent  heat  release  of  the 
clouds.  However,  the  adiabatic  heating  due  to  the  subsidence  outside  the  cloud 
cluster  was  found  to  contribute  significantly  to  the  evolution  of  Che  system. 
A  24-hour  integration  shows  that  the  model  is  capable  of  simulating  many 
important  features  such  as  Che  life  cycle,  the  intensity  of  circulation,  and 
rainfall  rates.  The  sensitivity  of  the  model  was  tested  by  including  the 
vertical  wind  shear,  deleting  Che  large-scale  forcing  at  the  middle  of  the 
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integration,  incorporating  the  cumulus  downdrafts  into  the  A-S  scheme,  and 
changing  the  mixed  layer  thermodynamic  properties.  Another  important  aspect  is 
that,  through  this  simulation  study,  we  were  able  to  observe  the  intrinsic 
nature  of  the  A-S  scheme  as  the  convective  system  evolved.  This  is  something 
which  those  previous  studies  did  not  emphasize.  For  comparison  purposes,  a  run 
with  the  use  of  Kuo's  1974  scheme  was  made.  It  was  found  that  Kuo's  scheme 
seems  unable  to  simulate  any  prominent  features.  We  believe  this  is  mainly  due 
to  the  fact  that  only  the  highest  clouds  are  allowed  to  develop  in  Kuo's 
scheme,  thereby  leading  to  the  sudden  Introduction  of  a  large  amount  of  heating 
to  the  upper  atmosphere.  Combined  with  the  cooling  effect  due  to  the  large- 
scale  lifting,  the  thermodynamic  structure  becomes  unfavorable  for  the 
development  of  the  cloud  cluster.  In  passing,  it  should  be  mentioned  that  the 
vertical  transport  of  horizontal  momentum  due  to  cumulus  clouds  is  not  Included 
in  the  model. 

The  topics  covered  in  this  report  range  from  the  investigation  of  the  A-S 
scheme  itself  to  the  applications  of  the  scheme  to  semi-prognostic  studies  and 
a  fully  prognostic  simulation.  Chapter  2  reviews  the  A-S  scheme  and  several 
related  studies.  Chapter  3  presents  the  results  of  the  semi-prognostic 
studies,  in  which  the  computational  algorithm  and  the  procedure  of  modifying 
the  A-S  scheme  by  including  the  downdraft  effects  are  also  discussed.  Chapter 
4  describes  the  prognostic  model  for  the  simulation  of  a  tropical  mesoscale 
convective  system.  Chapter  5  presents  the  model  results  and  sensitivity  tests. 
Chapter  6  gives  the  conclusion  of  this  report. 
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CHAPTER  II 


REVIEW  OF  THE  ARAKAWA-SCHUBERT  CUMULUS  PARAMETERIZATION 


In  this  chapter,  we  briefly  present  a  review  of  the  A-S  scheme  and  several 
previous  studies  in  which  the  A-S  scheme  was  applied.  The  full  details  of  the 
scheme  have  been  discussed  by  Arakawa  and  Schubert  (1974)  and  Lord  (1978).  The 
interested  reader  is  referred  to  those  original  articles.  The  large-scale 
budget  equations  are  reviewed  in  Section  2.1.  The  quasi-equilibrium  assumption 
of  the  cloud  work  function  Is  reviewed  in  Section  2.2 

2.1  The  large-scale  budget  equations; 


To  establish  the  relationship  between  the  collective  properties  of  the 
clouds  and  the  large-scale  variables,  the  heat  and  moisture  budgets  can  be 
expressed  in  z-coordinates  as; 


3  s  *  - 

P  +  pV*Vs 


+  pw  +  L(C  -  e) 


(2.1) 


p|i.p;.Vq.pw|I=(5-C)-i£|^ 


(2.2) 


where  the  variables  with  the  overbar  define  a  large-scale  system  and  the 
variables  with  the  prime  define  a  cloud-scale  system,  p  is  the  air  density, 
s  =  CpT  +  gz  is  the  static  energy,  q  is  the  mixing  ratio  of  water  vapor,  L  is 
the  latent  heat,  Qj^  is  the  radiative  heating  rate,  C  is  the  rate  of 


d  OW^  S  ^ 

condensation,  e  is  the  rate  of  evaporation  from  liquid  water,  and  — - -  and 

_  3z 

■'L-  respectively  represent  the  convergence  of  vertical  eddy  transports  of  s' 
3  z 

and  q'  with  clouds. 


Spectrally  dividing  the  cumulus  ensemble  into  sub-ensembles  according  to 
their  fractional  entrainment  rate  X  (this  variable  will  be  defined  shortly), 
the  vertical  eddy  transports  of  s'  and  q'  are  written  as: 


p  w' s" 


P  w'q' 


XjjCz) 

0 

Xn(z) 


m(2,X )( s^(z ,X)  -  s(z))dX 
m(z,X)(qj.(z,X)  -  q(z))dX 


(2.3) 


(2.4) 


where  the  subscript  c  denotes  a  cloud  variable.  The  notations  introduced  in 
(2.3)  and  (2.4)  will  be  explained  shortly.  Representing  the  sub-ensemble  cloud 
by  a  one-dimensional  steady-state  entraining  cloud  model  (the  governing 
equations  for  the  cloud  variables  are  given  in  Appendix  A),  the  vertical 
convergence  of  eddy  transports  of  s'  and  q'  can  be  expressed  as: 


3ow'  s' 
3z 


LC  -  M 
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c  3z 


(2.5) 


3pw'  q' 
3z 


D(q  -  q*)  -  Mj. 
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(2.6) 


The  detailed  derivations  of  (2.3)  through  (2.6)  are  given  in  Arakawa  and 
Schubert  (1974).  Substituting  (2.5)  and  (2.6)  into  (2.1)  and  (2.2), 


respectively,  we  obtain: 
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0  +  pV»7i  +  pw  -  -  Dli  +  M  |i  +  Qj.  (2.7) 

3  (  3  z  3  z 

p  M  +  pV»7q  +  pw  *  D(q*  -  q  +  £)  +  M.  (2.8) 

9 1  9  z  9  z 

A 

where  the  evaporation  rate  e  in  (2.1)  and  (2.2)  has  been  replaced  by  D1  whose 
physical  meaning  will  be  given  shortly.  Multiplying  (2.8)  by  L  and  adding  to 
(2.7)  gives  the  budget  equation  for  moist  static  energy  (h): 

p  +  pV«7h  +  pw  ■  DL(q*  ~  *1)  +  Qj.  (2.9) 

3  C  9  z  9  z 

where  h  ■  CpT  gz  +  Lq  and  the  asterisk  denotes  the  value  at  saturation  of  the 
variable.  The  notations  Introduced  in  (2.3)  to  (2.9)  are  explained  as  follows; 
(i)  ,  the  total  cloud  mass  flux,  is  expressed  as: 

M^(z)  -  m(z,X)dX  (2.10) 

where 

(2.11) 

m(z,X)dX  -  I  Mj^(z) 

^iC  (X  ,X+dX) 

is  the  sub-ensemble  mass  flux  due  to  the  clouds  which  have  the  parameter  Xj^  in 
the  interval  (X,X+dX).  Xq(z),  introduced  in  (2.3),  (2.4)  and  (2.10), 

represents  the  fractional  entrainment  rate  of  the  sub-ensemble  which  detrains 
at  the  cloud  top  at  the  level  z.  It  is  customary  in  diagnostic  and  prognostic 
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studles  of  cumulus  ensembles  to  divide  the  area  averaged  vertical  mass  flux 
(M),  which  can  be  expressed  as  pw.into: 

M=M^+M  (2.12) 

where  M  is  the  mass  flux  averaged  over  the  cloud-free  area. 


Each  individual  cumulus  cloud  is  considered  to  have  a  characteristic 
fractional  entrainment  rate  X  which  uniquely  determines  the  cloud  depth. 
Assuming  that  X  is  a  constant  with  height  for  each  cloud,  we  define  X  as: 

-—i -  3ni(_2_^_)_  ^  ^  (2.13) 

m(2,X)  32 

It  should  be  noted  that  the  total  cumulus  population  is  described  by  a 
continuum  of  cloud  entrainment  rates.  The  procedure  to  determine  X  is  given  in 
Appendix  A.  Solving  (2.13),  we  have 

m(2,X)  =  mg(X)  exp  {X(2  -  23)}  (2.14) 

where  mg(X)  is  the  cloud  mass  flux  at  cloud  base  and  2g  is  the  cloud  base 
height.  Here,  we  assume  that  2g  coincides  with  the  condensation  level  2^. 
From  (2.14),  we  define: 

I  exp  {X(2-2g)}  ,  2g  <  2  <  Zjj(X) 

(2.15) 


0 


,  2  >  2n(X) 
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as  Che  normalized  cloud  mass  flux. 

(li)  0(z),  Che  decrainmenc  of  cloud  mass  which  is  assumed  Co  cake  place  only 
aC  Che  cloud  Cop,  is  expressed  as: 


D(z) 


-m(z,Xjj(z)) 


dXgCz) 

dz 


(2.16) 


where  Xq(z)  denoces  Che  encrainmenc  raCe  of  Che  cloud  sub-ensemble  which 
decrains  ac  level  z. 

(Hi)  i(z)  is  Che  deCrained  liquid  waCer  which  spreads  at  Che  cloud  Cop  of  each 
sub-ensemble  and  evaporaces  Co  cool  Che  environmenc.  Il(z)  is  deCermined  from 
Che  one-dimensional  cloud  model,  as  shown  in  Appendix  A. 


The  physical  incerprecacions  of  Che  cloud  Cerms  in  (2.7)  and  (2.8)  are 

described  as  follows: 

3  s  3  Q 

(i)  M  —  and  M  —  are  respecCively  Che  cloud  hearing  and  drying  effeccs 

^  3z  3  z 

through  Che  forms  of  "cumulus- induced  subsidence."  Looking  inCo  (2.5)  and 

(2.6),  it  is  seen  that  the  real  physical  meanings  of  M  —  and  M  —  are  in 

3  z  ^  3  z 

terms  of  condensation,  convergence  of  eddy  transports,  and  the  detrainment  of 

water  vapor.  The  last  process  is  related  Co  M  —  only. 

3  z 

A 

(ii)  Di  is  the  evaporation  rate  from  the  detrained  liquid  water  at  the  cloud 
cop.  Without  the  evaporation  from  falling  precipitation  associated  with 
convective  downdrafts,  Dl  is  identical  to  e  in  (2.1)  and  (2.2). 

(iii)  D(q*-q)  is  the  detrainment  of  saturated  cloud  air  into  the  environment. 


It  is  customary  to  describe  (2.7)  and  (2.8)  by; 


Ql  =  P(|7  +  V.7i  ^  w  ||) 


■DU  +  —  +  Qr 


(2.17) 


Q,  =  -L  -t-  V*Vq  +  w  ii)  =1  -LD(q*  -  q  4)  -  LM^, 

^  9  ^  d  z  9  z 


(2.18) 


where  Qj^  and  Q2  are  respectively  Interpreted  as  the  "apparent  heat  source"  and 
"apparent  moisture  sink"  (Yanai  al.*  1  1973).  Yanal  £t  a^.  (1973)  first 
proposed  the  concept  that  the  heating  and  drying  effects  due  to  a  cloud 
ensemble  in  a  unit  area,  which  is  large  enough  to  contain  the  ensemble  but 
small  enough  to  cover  only  a  fraction  of  the  large-scale  disturbance,  can  be 
Indirectly  computed  from  large-scale  observational  data.  Since  the  emergence 
of  this  idea,  the  and  1^2  profiles  have  appeared  widely  in  numerous 
literatures  in  discussion  of  the  interaction  between  the  cumulus  clouds  and 
their  environment. 


To  complete  the  large-scale  budgets,  the  governing  equations  for  the. 
sub-cloud  mixed  layer  of  variable  depth  are  given  by  Arakawa  and  Schubert 
(1974)  as: 


“M  ^  ^  <rsv>„l  -  (Qr)„ 


(2.19) 


(2.20) 
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Equation  (2.21)  describes  the  time  change  of  the  depth  of  mixed  layer  (zg), 
which  implies  that  the  cloud  base  is  at  the  top  of  the  mixed  layer.  Equations 
(2.19)  and  (2.20)  give: 


-M  ^  .  k  ^  *  (Q„)^ 


(2.22) 


where  subscript  M  denotes  the  mixed  layer  values,  subscript  B  denotes  the  cloud 
base  values,  and  subscript  0  denotes  the  surface  values.  ,  F^,  Fj^,  and  F^^ 
are  the  vertical  turbulent  eddy  fluxes  of  s,  q,  h,  and  ,  respectively.  The 
symbol  t  represents  the  jump  of  a  particular  property  across  the  top  of  the 
sub-cloud  mixed  layer.  The  parameter  k  is  an  empirical  constant  which  is  the 
ratio  between  the  turbulent  eddy  flux  at  the  surface  and  the  value  at  the  top 
of  the  mixed  layer.  The  quantity  Mg  is  the  total  cloud  mass  flux  at  the  cloud 
base  and  Wg  is  the  large-scale  vertical  velocity  at  the  cloud  base.  Fig.  2,1 
is  a  schematic  diagram  of  the  sub-cloud  layer  thermodynamic  variables.  Arakawa 
and  Schubert  proposed  that  the  cumulus  clouds  originated  from  the  mixed  layer 
with  the  thermodynamic  values  of  and  qj,j.  A  consequence  of  this  assumption 
is  that,  as  seen  from  (2.19)  and  (2.20),  the  cumulus  convection  does  not 
directly  change  S|,j  and  qj^.  However,  S{,j  and  qj^  can  be  indirectly  modified  by 
cumulus  convection  through  the  change  of  the  depth  of  the  mixed  layer  (see 
(2.21);  Ogura  and  Cho,  197A;  and  Albrecht  aJL*  >  1979).  An  extensive 
investigation  of  the  coupling  processes  between  the  planetary  boundary  layer 
(PBL)  parameterization  and  the  A-S  cumulus  parameterization  in  a  GCM  is  given 


in  Suarez  et  al.  (1983).  A  detailed  discussion  along  this  line  will  be 


presented  in  Chapter  4,  where  we  develop  a  suh-cloud  layer  parameterization  for 
a  mesoscale  model. 


The  analysis  up  to  this  point  can  be  referred  to  as  the  first  part  of  the 
A-S  scheme.  Once  the  large-scale  variables,  including  —  and  — ,  are  given 

at  at 

from  observation,  (2.17)  and  (2.18)  can  be  used  to  diagnostically  determine  the 
cloud  properties  such  as  mg(X),  D(z),  and  i(z).  Ogura  and  Cho  (1973)  first 
demonstrated  the  properties  of  the  cumulus  population  by  applying  the  A-S 
scheme  to  the  mean  summertime  conditions  over  the  Marshall  Islands.  Their 
results  indicated  a  bimodai  distribution  of  the  cloud  base  mass  flux  in  terms 
of  the  height  of  the  cloud  top:  one  group  of  clouds  penetrated  up  to  the 
200-300  mb  level  and  another  group  stayed  below  the  600  mb  level,  while 
relatively  few  clouds  lay  between  400  and  600  mb.  Later,  Cho  and  Ogura  (1974) 
applied  the  same  method  to  the  averaged  easterly  wave  disturbances  in  the 
western  Pacific  based  on  Reed  and  Pecker's  data  set  (Reed  and  Recker,  1971). 
Simodal  distributions  are  also  present  in  the  region  near  trough  axis  while 
only  shallow  clouds  dominate  in  the  ridge  area.  Another  finding  is  that  there 
appears  to  be  a  linear  correlation  between  the  cloud  base  mass  flux  of  deep 
clouds  and  large-scale,  low-level  mass  convergence.  It  suggests  that  high 


cloud  activities  are  mostly  controlled  by  the  low-level,  large-scale  lifting. 


Nitta  (1975)  diagnostically  determined  the  cloud  mass  flux  distributions 
and  cloud  detrainment  from  observational  data  over  the  tropical  Atlantic  Ocean 
during  the  80MEX  (Barbados  Oceanographic  and  Meteorological  Experiment)  Phase 
III  (22-30  June  1969).  He  found  that  only  shallow  clouds  were  dominant  and 
that  the  total  cloud  mass  flux  was  confined  below  the  inversion  during  the 


Fig.  2.1  A  schematic  diagram  of  the  sub-cloud  layer  i hermodynamic  variables, 
and  are  the  mixed  layer  dry  static  energy  and  mixing  ratio  and  are 

sumed  constant  with  height,  h^  is  given  by  hj^  =  +  Lq^^.  (  )g  indicates 

a  value  in  the  free  atmosphere  immediately  above  the  PBL  top  at  z  =''’Zg.  There 
is  an  infinitely  thin  transition  layer  between  Zg  and  Zg  with  the  jump 
conditions  Aq,  As,  and  Ah. 
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typical  trade  iaversion  period.  The  cooling  and  moistening  effects  due  to  the 
detrainment  of  liquid  water  and  water  vapor  appear  to  be  very  crucial  in  the 
large-scale  heat  and  moisture  budgets  over  the  tropical  trade  wind  region. 
Based  on  the  framework  of  the  A-S  scheme,  Johnson  (1976)  included 
convective-scale  downdrafts  in  a  diagnostic  model  to  determine  the  cloud 
population  and  their  associated  properties.  Basically,  he  found  that  cumulus 
downdrafts  contributed  significantly  to  the  total  convective  mass  transport  in 
the  lower  atmosphere.  Johnson  concluded  that  the  neglect  of  cumulus  downdrafts 
and  their  associated  rainfall  evaporation  led  to  the  diagnosis  of  an 
extensively  large  population  of  shallow  clouds  in  highly  convective  situations. 
Johnson's  downdraft  model  will  be  discussed  in  detail  in  the  next  chapter, 
where  we  incorporate  his  model  into  the  A-S  scheme  in  a  prognostic  study. 
Later,  Johnson  (1980)  modified  his  diagnostic  model  by  including  mesoscale 
downdrafts  to  study  anvil-type  mesoscale  convective  systems.  An  entraining 
plume  model  with  a  spectrum  distribution  of  clouds  for  convective-scale 
elements  was  still  employed  as  before,  while  mesoscale  downdrafts  were 
formulated  through  a  bulk  cloud  model.  He  concluded  that  the  mesoscale 
downdrafts  are  equally  important  in  contributing  to  the  budgets  of  heat  and 
moisture . 


Undoubtedly,  using  the  A-S  scheme  to  diagnostically  determine  cloud 
properties  under  various  topical  weather  conditions  has  contributed  to  an 
understanding  of  the  interaction  between  cumulus  convection  and  large-scale 
disturbances.  Nevertheless,  the  ultimate  purpose  of  the  A-S  scheme  is  a 
prognostic  one,  and  the  tendency  terras  —  and  —  are  the  unknown  variables  to 

at  at 

be  determined  from  the  large-scale  variables.  This  problem  will  be  discussed 
in  the  next  section.  In  passing,  it  is  worthwhile  to  mention  a  physically 
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unreallstlc  aspect  which  may  be  encountered  when  using  the  A-S  scheme  in  a 
diagnostic  study.  With  the  use  of  (2.10)  and  (2. 140,  it  can  be  recognized  that 
(2.17)  and  (2.18)  are  a  set  of  two  Voiterra  integral  equations  of  the  second 
kind  with  two  unknowns  ll(z)  and  mg(X).  They  can  be  solved  by  a  proper 
numerical  approach  (Nitta,  1975).  However,  the  solutions  of  l{z)  and  mg(X) 
cannot  be  guaranteed  to  be  positive  values,  as  shown  by  Cho  and  Ogura  (1974) 
and  Johnson  (1976,  1978).  Although  the  magnitudes  of  the  negative  solutions  in 
their  results  are  relatively  small,  the  experience  gained  from  our  preliminary 
studies  shows  that  the  negative  magnitudes  of  Jl(z)  and  mg(X)  can  be  extensively 
large  in  some  cases  and  chat  these  situations  pose  difficulties  in  interpreting 
the  solutions. 

2.2  Cloud  work  function  and  quasi-equilibrium  assumption 

We  start  with  the  time  change  of  the  kinetic  energy  of  each  sub-ensemble: 

=  A(X)ra„(X)  -  D(X)  (2.23) 

dt  “ 

where  K(X)  and  D(X)  are  the  kinetic  energy  and  its  dissipation,  respectively. 
A(X)  is  called  the  cloud  work  function  and  is  a  measure  of  buoyancy.  The 
derivation  of  (2.23)  is  presented  in  Haltiner  and  Williams  (1979).  The  cloud 
work  function  A(X)  is  given  as: 

A(X)  =  /  - - -  n(z,X)(s  (z,X)  -  s(z)ldz 

CpT(z) 


(2.24) 


where  is  the  static  energy  of  the  sub-ensemble  X.  The  time  derivative 
A(X)  can  be  expressed  as: 


3A(X) 

at 


(2.25) 


fZp)(X)  ,Zn(^)  ah 

p(z)6(z)dz  +  p(z)b(2,X)Xn(z,X)  dz 


-  /  p(z)o(2)ri(z  ,X)  —  dz 


Substituting  the  expressions  of  and  from  (2.7)  and  (2.9)  into  (2.25), 
obtain: 


^1^  =  F|,,(X)  +  [K^(X,X')  +  Kjj(X,X')]m  (X')dX'  +  F  (X) 


(2.26) 


where  Fp(X)  =  f n (z ,X ) {Xb( z ,X )[ -w  ^  -  v*Vh  +  Q„] 


(2.27) 


F„(X) 


-  a(z)(-  V  —  -  vVs  + 

az 


-r—  /  p(z)6(z)dz 

at  •'  Zg 


Qr)}  p(z)dz 


(2.28) 


K^(x,X')  =  n(z,x)n(z,x'){xb(z,x)  —  -  0(2) 


(2.29) 
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Kij(X,X')  -  n(2j,,X)n(z5.X'){Xb(z5,X)  L(q*(25)  -  q(z5))  +  a(Z5)U(Z5) }  (2 .30) 
3 

where  — —  in  (2.28)  can  be  replaced  by  the  right  hand  side  of  (2.22)  or  another 
3 1 

expression  derived  from  a  different  approach  of  the  mixed  layer 
parameterization  (e.g.,  Suarez  et  al.,  1983  or  Johnson,  1976).  The  detailed 
derivations  of  (2.25)  and  (2.26)  and  the  actual  forms  of  a(z),  B(z),  and  b(z,X) 
are  presented  in  Appendix  B. 

The  differences  between  Equations  (2.26)  -  (2.30)  and  their  original 
expressions  in  Arakawa  and  Schubert  (1974)  are  due  to  the  following 
approximations  that  we  introduced: 

(i)  The  cloud  base  Zp  coincides  with  the  condensation  level  z_. 

H  c 

(ii)  The  effects  of  water  vapor  and  liquid  water  on  the  cloud  work  function 
are  neglected  in  (2.24). 

(lii)  The  top  of  the  mixed  layer  is  defined  by  a  fixed  z-coordinate  surface. 
Namely,  Zg  is  constant  with  time. 

The  first  two  approximations  are  adopted  from  Lord  ££  (1982)  for 

simplification  purposes.  The  last  approximation  may  not  be  valid  in  studying  a 
GCM  (Suarez  ^. ,  1983)  or  synoptic-scale  prediction  model.  However,  for  the 
serai-prognostic  studies  and  mesoscaie  simulations,  it  appears  to  be  acceptable 
(Lord,  1978,  1982;  Hack  £t  a^. ,  1983). 

The  kernel  K^(X,X')  in  (2.26)  serves  as  a  stabilization  factor  for  cloud 
type  X  due  to  the  heating  effect  of  cloud  type  X'  on  the  environment. 
Meanwhile,  F^(X)  and  (2.26)  serve  as  destabilization  factors  for  cloud 

type  X  due  to  large-scale  forcing.  Ku(X,X')  also  serves  as  a  destabilization 
factor  for  cloud  type  X  due  to  the  cooling  effect  through  the  evaporation  of 


detrained  liquid  water  and  the  moistening  effect  through  the  detrainment  of 
water  vapor  of  cloud  type  X'.  Arakawa  and  Schubert  (1974)  respectively  defined 
^c’  *^D  "cloud  layer  forcing,"  "mixed  layer  forcing,"  "vertical 
mass  flux  kernel,"  and  "detrainment  kernel."  However,  the  major  effects  in 
(2.26)  are  through  the  and  terms:  increases  the  cloud  work  function  of 
the  sub-ensemble  X  primarily  through  the  adiabatic  cooling  due  to  the  large- 
scale  upward  motion,  while  decreases  the  cloud  work  function  of  the  sub¬ 
ensemble  X  through  adiabatic  warming  of  the  environment  due  to  the  induced 
subsidence  of  sub-ensemble  X'.  Typically,  Fj^  and  are  respectively  smaller 
than  Fp  and  by  one  to  two  orders  of  magnitude.  However,  when  the  lower 
atmosphere  is  extremely  dry,  can  be  equally  as  important  as  due  to  the 
large  excess  mixing  ratio  (see  (2.30)).  Based  on  the  argument  of  the  quasi- 
equilibrium  assumption  (Arakawa  and  Schubert,  1974;  Lord  and  Arakawa,  1980), 
(2.26)  can  be  expressed  as: 

K(X  ,X")mg(X')dX'  +  F(X)  =  0,  when  mg(X)  >  0  (2.31a) 

K(X  ,X' )mg(X"  )dX'  +  F(X)  <  0,  when  mg(X)  =  0  (2.31b) 

Where  K(X,X')  =  1C^(X,X')  +  Kjj(X,X');  and  F(X)  =  F^(X)  +  Fj,j(X).  The  quasi¬ 
equilibrium  assumption,  which  serves  as  a  closure  hypothesis,  can  be  briefly 
Interpreted  as  follows:  if  the  adjustment  time  of  clouds  is  small  compared  to 
the  time  scale  of  the  large-scale  disturbance,  we  assume  that  the  cumulus 
ensemble  follows  a  sequence  of  quasi-equilibrium  states  while  the  large-scale 
forcing  changes;  and  the  condition  described  by  (2.31a)  must  be  valid  for  each 
positive  mg(X);  on  the  other  hand,  when  the  time  derivative  of  A(X)  becomes 
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negative,  a  sub-ensemble  of  such  X  would  not  exist,  thus  requiring  mg(X)  >  0  as 
shown  in  (2.31b). 

(2.31a)  is  a  Fedholm  integral  equation  of  the  first  kind  for  the  unknown 
mg(X)  subject  to  the  condition  of  Inequality  by  (2.31b).  Solving  (2.31)  is  not 
a  trivial  matter.  A  "simplex  over- adjustment  method"  proposed  by  Silva-Dias 
and  Schubert  (1977)  is  found  to  be  very  efficient,  which  is  based  on  a  linear 
programming  algorithm  (Dantzig,  1963;  Gass,  1975).  By  introducing  some  "slack" 
variables  into  (2.31a),  this  method  is  able  to  provide  positive  solutions  for 
mg(X)  and  meet  the  requirement  of  (2.3lb).  Generally  speaking,  the  linear 
programming  method  restates  the  quasi-equilibrium  hypothesis  as  an  optimization 
problem.  A  detailed  description  of  how  to  use  the  linear  programming  algorithm 
to  solve  (2.31)  is  given  by  Lord  (1978)  and  Hack  ad.  (1984). 

Once  mg(X)  is  determined,  M^(z)  can  be  readily  obtained  from  (2.14)  and 

(2.10)  and  D(z)  can  be  obtained  from  (2.16).  Substituting  M^,  D,  and  t  into 

(2.7)  and  (2.8)  with  the  given  large-scale  variables,  —  and  —  can  be 

d  t  9 1 

determined.  In  the  mean  time,  rainfall  rates  can  be  calculated  from  (A. 6). 

A  series  of  tests  of  the  A-S  scheme  have  been  conducted  by  Lord  (1978, 
1982),  Lord  and  Arakawa  (1980),  and  Lord  et  al.  (1982).  First,  Lord  and 
Arakawa  (1980)  assembled  several  data  sets  over  various  areas  in  the  tropics 
and  showed  the  evidence  that  the  cloud  work  function  is  a  quasi-uni versal 
function  of  the  detrainment  level  (i.e.,  the  height  of  cloud  top).  Later,  Lord 
(1982)  applied  the  semi-prognostic  approach  to  the  observational  data  for  the 
period  1-18  September  1974  over  the  GATE  area  to  evaluate  the  A-S  scheme.  The 
predicted  vertical  profile  of  the  time-averaged  heating  effects  of  clouds 


'  v.v.v.' 


agrees  well  with  Che  observed  one  (Fig.  2.2a),  while  Che  drying  effecCs  of 
clouds  show  less  saCisfacCory  agreemenc  (Fig.  2.2b).  Also,  Che  predicCed 
rainfall  races  agree  quice  well  wich  Che  observacions . 


The  algorichm  which  Lord  (1982)  used  in  his  semi-prognosCic  scudy  can  be 
ouclined  in  che  following  seeps: 

(i)  From  che  large-scale  variables  T,  q,  w,  u,  v,  and  Q^,  cloud  parameCers  X, 
n(z,X),  h^,  s^,  and  q^  are  compuCed  for  each  sub-ensemble  whose  cloud  Cop 
heighc  is  specified  according  Co  Che  vercical  grid  arrangemenc  of  che  model 
aCmosphere. 

(ii)  The  cloud  work  funecion  (2.24)  is  compuCed  for  each  sub-ensemble. 

(iii)  The  large-scale  forcing  in  (2.31)  for  che  ich  sub-ensemble  is  compuCed 
by: 


F(i) 


A^(i)  -  A(i) 
AC' 


(2.32) 


where  A(i)  is  calculaCed  from  Che  daCa  ac  a  given  observaCion  cime  c  and  A'(c) 
is  che  cloud  work  funeCion  compuCed  from  bogus  large-scale  cemperaCure  and 
moisCure  fields  which  would  resulC  when  Che  large-scale  horizoncal  and  vercical 
adveccion  operaCed  for  a  shore  period  of  Cime  AC'.  NoCe  ChaC  Che  fracCional 
enCrainmenC  race  for  che  iCh  sub-ensemble  in  A'(i)  is  differenc  from  ChaC  in 
A(i)  because  Che  CemperaCure  and  moisCure  profiles  have  been  modified. 

(iv)  The  kernel  in  (2.31)  is  compuCed  by: 
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Fig.  2.2.  (a)  the  time  averaged  vertical  distributions  of  the  apparent  heat 
source  due  to  clouds  (Qi“Qd)  The  time  averaged  vertical  distributions 
of  the  apparent  moisture  sink  (Qo)*  Boch  are  in  the  units  of  C  day  .  Ihe 
predicted  values  are  indicated  by  solid  lines  and  the  observed  values  are 
indicated  by  dashed  lines.  (Adapted  from  Lord,  19B2.) 


where  A  (1)  is  computed  after  the  large-scale  variables  are  modified  by  the  jth 
sub-ensemble  through  its  cloud  effects  in  the  time  interval  At".  Again,  the  X 
in  A"(i)  is  different  from  the  original  one  in  A(i). 

(v)  Fredholm's  integral  equation  (2.31)  is  now  solved  for  mg(X)  by  the 
"simplex  over-adjustment  method". 

(vi)  —  and  rainfall  rate  are  computed  from  (2.7),  (2.8),  and  (A. 6), 
respectively. 


It  is  seen  that  Lord  used  finite  difference  approximations  in  steps  (iii) 
and  (iv)  to  calculate  forcing  and  kernel  terms  instead  of  using  their  explicit 
forms  stated  in  (2.27)  to  (2.30).  Further  discussion  of  Lord's  algorithm  is 
given  in  Chapter  3.  Recently,  Hack  et  al.  (1984)  modified  the  original  A-S 
scheme  into  a  flux  form  for  computational  convenience.  They  subsequently 
applied  the  scheme  to  a  semi-prognostic  study  for  the  Marshall  Islands  data  and 
a  prognostic  study  of  the  development  of  an  axisymmetric  hurricane.  However, 
they  focused  their  attention  primarily  on  the  numerical  aspects  of  the  A-S 
scheme  without  much  discussion  of  the  model  results.  The  only  conclusion  they 
presented  about  the  physical  aspects  of  their  study  is  that  the  mass  flux 
associated  with  deep  clouds  is  controlled  not  only  by  vertical  motion  at  the 
top  of  the  mixed  layer  but  also  by  large-scale  vertical  motion  in  the  middle 


and  upper  troposphere. 


CHAPTER  III 


THE  SEMI-PROGNOSTIC  STUDIES 

In  this  chapter,  we  present  the  results  of  semi-prognosticaily  applying 
the  A-S  scheme  to  two  different  data  sets:  one  for  an  intense  cloud  band 
observed  in  GATE,  and  the  other  for  wave  composites  of  GATE  Phase  II  and  III. 
Section  3.1  presents  the  motivation  of  this  study.  Section  3.2  gives  the 
computational  algorithm.  Section  3.3  presents  the  results.  Section  3.4 
describes  an  experiment  which  incorporates  the  convective-scale  downdraft 
effects  into  the  A-S  scheme  for  prognostic  purposes. 

3.1  Introduction: 

Although  the  series  of  studies  by  Lord  (1978);  Lord  and  Arakawa  (1980), 
Lord  (1982),  and  Lord  £t  £l.  (1982)  have  been  considered  the  most  extensive 
approach  to  test  the  A-S  scheme,  we  feel  that,  as  far  as  measuring  the  success 
of  the  A-S  scheme  through  a  semi-prognostic  approach  is  concerned,  simply 
presenting  the  result  averaged  over  the  period  of  18  days  (1-18  September  1974 
of  GATE  Phase  III),  as  shown  in  Fig.  2.2,  may  not  be  an  accurate  assessment  of 
the  scheme.  In  order  to  examine  the  response  of  cumulus  clouds  to  the  large- 
scale  forcing  through  the  application  of  the  A-S  scheme  to  an  "individual" 
tropical  convective  event,  a  semi-prognostic  approach  using  the  data  for  an 
Intense  cloud  band  observed  from  the  GATE  A/B  scale  ship  array  on  12  August 
1974  was  undertaken.  Furthermore,  in  order  to  examine  the  general  behavior  of 
the  A-S  scheme  in  the  mean  weather  conditions  over  the  tropics,  we  applied  the 
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scheme  to  the  wave  composites.  The  composite  data  of  GATE  Phase  II  and  III 
were  used  for  this  purpose. 

To  date,  all  the  studies  on  the  A-S  scheme  have  placed  the  emphasis  on  the 
net  heating  and  drying  effects  of  cumulus  clouds,  neglecting  their  associated 
individual  components  ,  such  as  condensation,  evaporation,  and  vertical 
transports  of  heat  and  moisture,  which,  in  fact,  can  be  easily  deduced  from  the 
A-S  scheme.  Soong  and  Tao  (1980),  Tao  (1983),  and  Soong  and  Tao  (1984)  applied 
the  "cloud  ensemble  model"  developed  by  Soong  and  Ogura  (1980)  to  the  GATE 
cloudband  case.  Their  results  indicate  that  statistical  properties  of  model 
clouds  that  develop  under  a  given  large-scale  forcing  agree  well  with  the 
observations  in  many  aspects.  We  shall  compare  not  only  the  predicted  cloud 
heating  and  drying  effects  with  the  observations  but  also  the  predicted  rates 
of  condensation,  evaporation,  and  vertical  transports  of  heat  and  moisture  with 
those  derived  from  Soong  and  Tao's  cloud  modeling  study.  This  comparison 
enables  us  to  further  illustrate  the  physical  details  of  the  A-S  scheme. 

It  has  become  increasingly  evident  that  in  regions  of  active  convection 
the  convective-scale  downdrafts  and  their  associated  evaporation  from 
precipitation  are  important  not  only  during  the  life  cycle  of  individual 
cumulus  clouds  but  also  in  contributing  to  the  budgets  of  large-scale 
disturbances.  On  the  other  hand,  the  comparative  study  mentioned  before 
reveals  that  the  A-S  scheme  underestimates  both  condensation  and  evaporation 
rates  substantially  when  compared  with  the  results  form  the  cloud  ensemble 
model.  We  believe  that  the  inclusion  of  the  convective-scale  downdraft  effects 
may  alleviate  this  deficiency.  As  mentioned  in  Chapter  2,  Johnson  (1976) 
developed  a  method  to  diagnostically  estimate  the  contribution  to  cloud 
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properCles  by  convective-scale  downdrafts.  Based  on  Johnson's  model,  a 
formulation  which  incorporates  the  downdraft  effects  Into  the  quasi-equilibrium 
assumption  is  developed  and  applied  to  the  semi-prognostic  study.  The  result 
shows  that  a  substantial  improvement  of  predicted  condensation  and  evaporation 
is  obtained. 

3.2  The  algorithm: 

Because  the  cumulus  parameterization  is  incorporated  in  a  vertically 
discrete  model  atmosphere,  the  use  of  the  fractional  entrainment  rate  X  as  the 
spectral  parameter  in  the  original  theory  of  the  A-S  scheme  may  introduce 
computational  difficulties.  Namely,  if  we  specify  X  in  time  as  an  independent 
variable,  then  the  detrainment  level  Zj)(X)  will  be  a  function  of  time,  which  is 
not  compatible  with  the  fixed  vertical  coordinates.  Thus,  in  the  applications 
of  the  A-S  scheme  we  may  determine  the  cloud  spectrum  distributions  with  the 

A  A 

use  of  cloud  top  height,  denoted  by  z,  instead  of  using  X.  z  is  associated 
with  the  pre-existing  vertical  grids  and  X,  which  corresponds  to  each  z,  is 
determined  by  an  iterative  method  through  the  large-scale  thermodynamic 
structure  with  the  cloud  top  condition  of  zero  buoyancy  (see  Appendix  A). 
Hence,  X  is  no  longer  an  independent  variable.  In  this  way,  we  can  rewrite 
(2.31)  as: 

A 

l^max  K(z  ,z' )mg(z' )dz'  +  F(z)  =•  0  ,  when  mg(z)  >0  (3.1a) 

B 

j^max  K(z,V)n,g(z')d2'  +  F(z)  <  0  ,  when  mg(z)  =  0  (3.1b) 

B 

where  the  expressions  for  k(z,  ) ,  and  F(z)  remain  the  same  as  shown  in  (2.27) 


to  (2.30).  Note  that  the  upper  and  lower  limits  of  the  integral  equations  have 

been  changed  according  to  the  new  spectral  parameter.  In  transforming  (2.31) 

to  (3.1),  all  the  terms  associated  with  the  time  rate  of  change  of  X  f — )  were 

at-^ 

neglected.  As  indicated  by  Hack  et  al.  (1984),  where  the  pressure  at  the 
cloud  top  (p)  was  used  to  represent  the  cloud  spectrum,  this  change  in  spectral 
parameter  somewhat  alters  the  quantitative  details  of  quasi-equilibrium  as 

ax 

originally  formulated  because  of  the  negligence  of  —  in  deriving  (3.1a)  and 

a  t 

(3.1b).  Nonetheless,  the  qualitative  aspects  should  remain  unchanged,  since 
the  modifications  of  the  kernel  and  forcing  terms  due  to  the  time  rate  of 
change  of  X  is  small  (Schubert,  1973).  A  complete  discussion  of  the  use  of  the 
detrainment  level  (z,  or  p)  as  the  cloud  spectral  parameter  in  the  A-S  scheme 
is  given  in  Hack  ^  a^.  (1984),  where  interested  readers  should  be  referred 
to.  From  now  on,  all  the  Xs  in  the  equations  related  to  the  A-S  scheme  are 
replaced  by  the  new  spectral  parameter  z. 

It  was  noticed  in  Section  2.2  that  Lord  used  the  finite  difference 
approximations  (2.32)  and  (2.33)  to  calculate  the  forcing  and  kernel  terms 
instead  of  using  the  explicit  forms  as  stated  in  (2.27)  to  (2.30).  The 
advantage  of  Lord's  approach  is  that  the  time  rate  of  change  of  X  can  be 
numerically  determined  with  two  bogus  time  intervals  At'  and  At"  and  a  bogus 
mass  flux  mg(j).  In  principle.  At'  in  (2.32)  and  mg(j)At"  in  (2.33)  can  be 
arbitrarily  chosen  (Lord,  1978).  However,  it  was  indicated  by  Lord  that  there 
are  practical  upper  and  lower  bounds  for  the  choices.  So  (1982)  applied  the 
A-S  scheme  using  Lord's  algorithm  to  the  mid-latitude  Baiu  situation.  He  found 
that  A"  -  A  in  (2.33)  is  not  linear  with  the  choice  of  mg(j)At"  so  that  the 
kernel  and,  as  a  consequence,  the  cloud  base  mass  flux  (nig)  depend  upon  the 
choice  of  mg(J)At".  Figs.  3.1  and  3.2,  adapted  from  So  (1982),  indicate  that 
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Fi2.  3.1.  The  time  series  of  .cloud  base  mass  fluxes  in  the  units  of  mb  r- 
for  (a)  mb(j)at"  =  50^Kg  m  ^  and  (b)  mb(j)&t"  =  85  Kg  m  .  The  ordinates  are 

the  cloud  top  pressure  p  (mb).  (Adapted  from  So,  1982.) 


there  are  substantial  differences  both  in  nig  and  the  precipitation  rates  when 
two  different  mg(J)dt"  are  used.  Our  preliminary  examination  also  reveals  that 
the  different  choices  of  At'  and  mg(j)At"  result  in  different  heating  and 
drying  profiles.  Although  the  differences  are  not  as  dramatic  as  So  showed,  we 
feel  that  a  more  acceptable  algorithm  is  required. 

To  avoid  the  undesirable  sensitivities  described  above,  we  applied  the  A-S 
scheme  by  directly  calculating  Che  forcing  and  kernels  from  (2.26)  to  (2.29), 
but  with  the  change  of  spectral  parameter  from  k  to  z.  As  mentioned  in  Chapter 

A  A  A 

2,  once  F(2)  and  K(z,z')  in  (3.1)  are  determined,  mg  can  be  obtained  by  solving 
an  optimization  problem.  The  explicit  forms  of  this  problem  can  be  Illustrated 
as  follows: 

Numerically  using  the  summation  to  represent  the  integral,  (3.1a)  can  be 
expanded  as: 

a(l,l)X(l)  +  a(l,2)X(2)  +  .  +  a(l,I^3x)  X(I^^^)  =  b(l) 

a(2,l)X(l)  +  a(2,2)X(2)  +  .  +  a(2,I^ax)  X(I^^^)  =  b(2) 

...  ...  (3.2) 


where  I^ax  denotes  Che  number  of  total  sub-ensembles,  a(n,m)  =  -K(z,z'), 
b(n)  =■  F(z)  and  X(m)  =  mo(z')dz'  for  n,m  =  1,  2,  .... 

o  UiaX 

Introducing  the  other  "slack"  variables  expressed  as  X(I  +1),  X( I  +2) 

HiaX  ni3X  nidx 

.  into  (3.2)  with  the  negative  signs  in  front  of  them,  we  have; 
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a(l,l)X(l)  +  a(l,2)X(2)  +  . +  a(  1 , I^ax^  "  +  !>  =  b(l) 

a(.M)X(l)  +  a(2,2)X(2)  +  . +  a(2 , I„ax)X( -  XCl^^^  +  2)  =  b(2) 

(3.3) 

•  •  •  •  •  • 

^ . W«W  - 

Now,  we  can  rewrite  (3.3)  into  the  form  of  the  "Linear  Programming  Problem" 
(also  known  as  an  optimization  problem)  as  follows:  Solve  the  set  of  linear 
equat ions : 

2x1 

„max 

y  a(n,m)  X(m)  =  b(n),  for  n  =  I,  2,....  (3.4) 

^  iTlaX 

ni=l 

subject  to  the  constraints  of 

X(m)  >  0  for  m  =  I,  2,  . 2x1^3,^  (3.5) 

and  a  minimized  linear  objective  function  Z,  defined  by: 

2x1 

„max 

Z  =  I  c(m)X(m)  for  m  =  1,  2 . .  2x1^^^  (3.6) 

m=l 


where 


c(m)  =  0 ,  for  m  =  1  ,  2 


V 


(3.7a) 
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Extenslve  discussions  of  this  linear  programming  problem  are  given  in  Lord 
(1978)  and  Hack  al .  (198A).  Our  biggest  concern  has  been  to  obtain  a 
software  package  which  can  provide  the  most  accurate  solutions.  It  was  found 
that  a  subroutine  named  ZX4LP,  based  on  the  work  by  Hanson  and  Wisniewski 
(1979),  from  the  IMSL  (International  Mathematic  and  Statistical  Library) 
computing  package  is  very  efficient  and  accurate  enough  for  our  purposes. 


Our  computational  algorithm,  shown  in  Fig.  3.3,  is  generally  the  same  as 
Lord's  described  in  Section  2.2  except  that  the  calculations  for  the  forcing 
and  kernel  terms  are  different.  In  this  study,  we  use  z-coordinates  and  divide 
the  atmosphere  vertically  into  16  layers  with  the  equal  depth  of  1  km.  As 
stated  in  Section  2.2,  is  in  general  much  smaller  than  F^.  (see  Eqs.  (2.26), 
(2.27),  and  (2.28)).  Consequently,  an  additional  assumption  that  hj^^  is 
independent  of  time  (i.e.,  *  0)  is  used  in  this  semi-prognostic  study  for 
simplicity.  We  followed  Lord's  (1982)  discrete  analogues  of  the  A-S  scheme. 
Fig.  3.4  shows  the  vertical  structure  of  the  large-scale  variables.  The 
schematic  diagrams  for  the  arrangement  of  the  cloud  variables  and  the 
distribution  of  the  cloud  spectrum  are  respectively  shown  in  Figs.  A.l  and  A. 2. 


3.3  The  results: 


3.3.1  The  tropical  cloud  band  study: 

A  complete  data  analysis  for  the  intense  cloud  band  that  developed  on  12 
August  1974  over  the  GATE  A/B  scale  ship  array  is  given  by  Ogura  et  al . 
(1979).  The  radar  photographs  shown  in  Fig.  3.5  indicate  that  this  system 
started  developing  around  00  GMT  and  reached  its  peak  activity  around  15  GMT. 
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Fig.  3.3.  The  schematic  diagram  of  the  computational  procedure  of  the  A-S 
scheme  in  the  semi-prognostic  study. 
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Fig.  3.4.  The  vertical  structure  of  the  large-scale  model.  Pressure  (p)  is 
provided  at  both  integer  and  half-integer  levels,  temperature  (T)  and  mixing 
ratio  (q)  are  provided  at  integer  levels,  and  vertical  velocity  (w)  is  given  at 
half-integer  levels.  The  thermodynamic  variables  at  half-integer  levels 
T(k+l/2)  and  q(k+l/2)  are  the  interpolated  values  from  those  at  k+1  and  k. 


trlous  convective  stages  at  (a)  00  GMT,  (b) 


KlK.  Radar  photographs  at  vai 

06  OMT.  (c)  12  GMT,  (d)  18  GMT,  and  (e)  24  GMT  on  12  August  1974. 
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Their  results  indicated  that  the  formation  of  this  convective  system  was 
preceded  by  mesoscale  low-level  convergence.  Later,  Soong  and  Tao  (1980, 
1984),  and  Tao  (1983)  used  the  cloud  ensemble  model  (Soong  and  Ogura,  1980)  to 
examine  the  response  of  the  cloud  band  to  the  mesoscale  lifting  processes 
described  by  Ogura  a^.  (1979). 

We  have  applied  the  A-S  scheme  to  this  tropical  convective  system  using  a 
semi-prognostic  approach.  Three-hourly  rawinsonde  observations  and 
Icinematically-computed  vertical  velocities  were  available  for  this  event.  The 
radiative  cooling  rate  (Q^)  used  in  this  study  (shown  in  Fig.  3.6)  is  the 
climatological  profile  analyzed  by  Dopplick  (1970).  Figs.  3.7  and  3.8  show  the 
results  of  predicted  and  Q2  respectively.  By  "predicted"  -  Qj^  and 
Q2,  we  mean  those  terms  associated  with  cloud  effects  in  the  last  members  of 
(2.17)  and  (2.18).  For  comparison  purposes,  the  profiles  of  and  Q2 
estimated  from  the  observed  large-scale  heat  and  moisture  budgets  are  also 
shown.  Agreement  between  the  predictions  and  observations  is  generally  good, 
particularly  for  the  heating  profiles.  We  note  that,  in  contrast  to  Fig.  2.2 
which  represents  the  results  averaged  over  the  period  of  18  days,  the 
predictions  shown  in  Figs.  3.7  and  3.8  are  for  the  observations  of  3-hour 
Intervals . 

Fig.  3.9  shows  the  total  mass  flux  (M)  (cf.  (2.12))  as  a  function  of  time 
and  height.  The  evolution  of  the  system  in  terms  of  M  is  evident.  It  is  seen 
that  a  low-level  maximum  is  present  during  the  period  of  03  to  06  GMT,  whereas 
a  high-level  maximum  is  present  during  09  to  18  GMT.  Fig.  3.10  shows  the  cloud 
mass  flux  (M^)  as  a  function  of  time  and  height.  It  is  observed  that  the 


general  feature  of  the  evolution  of  M  shown  in  Fig.  3.10  is  very  similar  to 
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Fig,  3.7.  The  vertical  distributions  of  the  apparent  heat  source  due 
(Qi-Qr)  at  (a)  03  z,  (b)  06  z,  (c)  09  z,  (d)  12  z,  (e)  15  z,  (f)  18 
21  z.  The  predicted  values  are  indicated  by  solid  lines,  and  the 
values  are  indicated  by  dashed  values. 
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that  in  Fig.  3.9  both  qualitatively  and  quantitatively,  implying  that  M  (cf. 
(2.12))  is  negligibly  small.  This  aspect  is  supported  by  the  cloud  model 
result  simulated  by  Soong  and  Tao  (1980,  1984),  which  will  be  discussed 
extensively  in  Section  3.3.3.  Fig.  3.11  shows  the  cloud  base  mass  flux  (mg)  as 
a  function  of  cloud  top  height  and  time.  The  results  indicate  that  only  deep 
clouds  with  cloud  top  height  above  8  km  (350  mb)  prevailed  over  the  mature  and 
decaying  periods,  with  the  maximum  around  12  to  15  GMT,  while  both  shallow  and 
deep  clouds  were  present  for  the  developing  period  from  03  to  09  GMT,  making  a 
bimodal  distribution  for  this  period.  This  implies  that  the  shallow  clouds  are 
associated  with  the  low-level  maximum  of  the  total  mass  flux  (M). 

Figs.  3.12  and  3.13  show  the  vertical  profiles  of  the  contributions  of 
individual  components  in  the  last  members  of  (2,17)  and  (2.18)  to  the  predicted 
(Qj  -  Qg)  and  Q2  at  09  GMT.  As  expected,  the  effects  of  cumulus- induced 
subsidence  are  the  dominant  factors.  However,  at  this  particular  time,  it  is 
noticed  that  the  low-level  detrainment  of  water  vapor  is  also  significant  in 
contributing  to  the  moisture  budget.  In  general,  this  effect  is  significant 
when  shallow  clouds  are  abundant  and  the  atmosphere  is  relatively  dry  at  the 
lower  levels.  More  discussion  will  be  presented  in  Section  3.3.3  where  we 
compute  the  cloud  properties  in  terms  of  condensation,  evaporation,  and 
vertical  eddy  transports  of  heat  and  moisture  and  compare  them  against  the 
results  from  the  cloud  ensemble  model. 

3.3.2  Easterly  wave: 

Chen  and  Ogura  (1982)  have  prepared  a  composite  of  eight  tropical  wave 
disturbances  that  passed  over  the  GATE  A/B  ship  array.  The  filtered  700  mb 
meridional  wind  was  analyzed  and  one  of  the  eight  wave  categories  was  assigned 
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Fig.  3.12.  The  vertical  distributions  of  the  contributions  due  to  cumulus- 

3  s 

induced  subsidence  and  evaporation  of  the  detrained  liquid  water  (-Di) 

to  the  predicted  net  heating  rate  indicated  by  "p".  The  net  heating  rate 
estimated  from  large-scale  observations  Is  indicated  by  "0". 
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Fig.  3.13.  The  vertical  distributions  of  the  contributions  due  to  cumulus- 
induced  subsidence  (M  ^),  detrainment  of  water  vapor  (D(q-q)  and  evaporation 

A 

of  detrained  liquid  water  (Di)  to  the  predicted  cloud  drying  effect  indicated 


by  "p".  The  net  cloud  drying  effect  estimated  from  large-scale  observation  is 
indicated  by  "0". 


to  each  observation  time  at  the  center  of  the  array  (i.e.,  23.5  W,  8.5  N)  in 


three  phases  of  GATE.  The  rawinsonde  and  surface  data  from  the  A/B-  and 
B-array  ships  were  then  composited  according  to  the  wave  categories.  Wave 
phases  or  categories  have  been  numbered  1  to  8 ,  with  category  2  centered  on  the 
region  of  maximum  northerly  wind,  category  4  on  the  trough,  category  6  on  the 
region  of  maximum  southerly  wind  and  category  8  on  the  ridge.  Categories  1,  3, 
5  and  7  occupy  intermediate  positions. 

In  order  to  examine  the  general  behavior  of  the  A-S  scheme  in  various 
weather  conditions  in  the  tropics,  composited  easterly  wave  data  in  Phase  II 
acd  III  were  used  in  this  semi-prognostic  study.  Phase  I  was  not  included, 
because  it  is  not  very  reliable  (Chen,  private  communication).  The  total  mass 
flux  (M)  is  shown  in  Fig.  3.14.  Significant  features  of  M  in  Phase  II  are  the 
middle  Level  maxima  around  categories  2  and  3,  low  level  maxima  around 
categories  5  and  6, and  relatively  strong  downward  motion  around  categories  1 
and  8.  On  the  other  hand,  in  Phase  III,  only  low  level  maxima  are  present  in 
categories  3,  4,  and  5,  whereas  very  weak  low-level  upward  motion  is  present  in 
categories  7  and  8.  It  is  also  seen  that  very  weak  downward  motion  is  present 
around  middle  and  upper  levels  in  categories  7  and  8  in  Phase  III.  Figs.  3.15 
and  3.16  present  the  observed  Qj  -  Qj^  and  Q2 ,  respectively,  whose  general 
features  are  similar  to  those  of  M.  The  predicted  Qj^  -  and  Q2  are 
respectively  shown  in  Figs  3.17  and  3.18.  Comparing  Figs.  3.17  and  3.18  with 
Figs.  3.15  and  3.16,  respectively,  it  is  found  that  the  agreement  between  the 
prediction  and  observation  is  generally  good.  No  cloud  heating  and  drying 
effects  exist  in  categories  1  and  8  in  Phase  II  simply  because  of  the  presence 
of  strong  downward  motion  in  those  categories.  The  solutions  for  the  cloud 
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Fig.  3.14.  The  total  mass  fluxes  (M)  in  the  units  of  mb  hr”^  for  the  composite 
waves:  (a)  for  Phase  II  and  (b)  for  Phase  III  as  functions  of  wave  category 
and  height. 
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Fig.  3.15.  The  observed  -  Qj^  in  the  units  of  C  day~^  for  the  composite 
waves:  (a)  for  Phase  II  and  (b)  for  Phase  III  as  functions  of  wave  category 
and  height. 
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Fig*  3 •16*  The  observed  Q2  in  the  units  of  C  day  ^  for  the  composite  waves i 
(a)  for  Phase  II  and  (b)  for  Phase  III  as  functions  of  wave  category  and 


IT 


Fig.  3.18.  The  predicted  Q 
(a)  for  Phase  II  and 
height . 


in  the  units  of  C  day”^  for  the  composite  waves: 

as  functions  of  wave  category  and 


(b^  for  Phase  III 


M 


base  mass  flux  as  functions  of  cloud  top  height  and  category  are  shown  in 
Fig.  3.19.  The  bimodal  distributions  are  clearly  shown  in  categories  3,  4,  and 
3  for  Phase  II  and  in  categories  4  and  5  for  Phase  III,  where  strong  convective 
activities  are  present.  This  also  reveals  the  relationship  between  the  shallow 
clouds  and  the  low-level  maximum  of  M  as  described  in  the  preceding  section. 
On  the  other  hand,  only  shallow  clouds  dominate  around  categories  7  and  8  for 
Phase  III  because  upward  motion  is  confined  to  below  6  km  in  those  categories. 
The  total  cloud  mass  flux  is  shown  in  Fig.  3.20.  It  is  also  seen  that  the 
general  features  and  magnitudes  of  are  similar  to  those  of  M,  shown  in 
Fig.  3.14. 

We  have  presented  the  results  of  a  semi- prognostic  application  of  the  A-S 
scheme  to  the  various  data  sets  over  the  GATE  area.  This  semi-prognostic  study 
is  primarily  intended  to  demonstrate  the  tractability  and  potential  advantages 
of  the  quasi-equilibrium  assumption,  which  is  considered  the  core  of  the  A-S 
scheme.  We  are  encouraged  by  the  fact  that  this  scheme  not  only  physically 
permits  a  highly  interactive  behavior  between  the  large-scale  disturbances  and 
the  cumulus  ensemble,  but  also  gives  realistic  predictions  of  the  cloud  heating 
and  drying  effects. 

3.3.3  The  comparison  with  the  cloud  modeling  study 

Soong  and  Ogura  (1980)  developed  a  "cumulus  ensemble  model".  Basically, 
it  is  a  nonhydrostatic,  anelastic  moist  cloud  model  with  parameterized  cloud 
physics  and  turbulent  processes.  Using  this  model,  one  can  study  the 
statistical  feedback  effects  of  a  group  of  cumulus  clouds  on  the  large-scale 


environment  in  which  a  group  of  clouds  is  embedded.  Soong  and  Tao  (1980,  1984) 
and  Tao  (1983)  have  applied  this  model  to  the  cloud  band  described  previously. 
The  thermodynamic  feedback  effects  of  cumulus  clouds  on  the  mesoscale 
temperature  and  mixing  ratio  profiles  due  to  condensation,  evaporation,  and 
vertical  eddy  transports  of  heat  and  moisture  with  clouds  were  investigated. 
This  experiment  enables  us  to  gain  some  insight  into  the  underlying  physical 
processes  and  quantitative  aspects  of  the  cloud  properties  which  cannot  be 
directly  measured.  On  the  other  hand,  this  model  may  be  categorized  as  a 
"quasi-prognostic"  approach,  since  the  large-scale  advectlon  of  temperature  and 
moisture  in  the  model  is  prescribed  from  the  observational  data.  Thus,  the 
comparison  between  its  results  and  those  deduced  from  the  semi-prognostic  study 
described  in  Section  3.3.1  should  be  physically  sound. 

Adapted  from  Tao  (1983),  the  profiles  of  simulated  and  observed  cloud 
heating  effects  averaged  in  the  period  from  06  to  12  GMT,  which  was  the 
developing  stage  of  the  cloud  band,  are  shown  in  Fig,  3.21.  The  results  from 
the  semi-prognostic  study  at  09  GMT  has  been  shown  in  Fig.  3.7c.  Both  results 
are  encouraging  in  that  the  model  cloud  ensemble  reaches  a  quasi-equilibrium 
state  in  which  the  heating  and  drying  effects  of  clouds  are  almost  balanced  by 
the  cooling  and  moistening  effects  due  to  large-scale  advection  and  radiative 
cooling,  as  should  be  the  case  in  the  real  atmosphere.  One  of  the  important 
conclusions  derived  by  Soong  and  Tao  (1980)  and  Tao  (1983)  is  that  (the 
cloud  mass  flux)  is  very  close  to  M  (the  total  mass  flux)  as  shown  in 
Fig.  3.22.  As  already  shown  in  Figs.  3.9  and  3.10,  the  M^  derived  by  the  A-S 
scheme  is  also  very  close  to  M. 


In  the  A-S  scheme  the  cloud  effects  due  to  condensation  rate  (C)  and 


vertical  convergence  of  eddy  transports  of  static  energy  and  moisture  ( 
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and  are  implicitly  included  in  the  terms  of  cumulus-induced  subsidence 

d  z 


(M„—  and  ) .  This  can  be  verified  by  the  relationship  among  (2.1),  (2.5), 
C3s  ^3z  _  _ 


and  (2.7).  Meanwhile,  from  (2.3)  and  (2.4),  ^  and  can  be  expressed 
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m(z,z)(q  (z,z)  -  q(z))dz 
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Hence,  once  mD(z)  is  known. 


3pw' s" 


and 


3pw'q' 


can  be  obtained  from  (3.8)  and 
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(3.9),  respectively.  Then,  we  can  calculate  C  either  from  (2.5)  or  (2.6).  On 

the  other  hand,  as  indicated  in  Secton  2.1,  the  evaporation  rate  (e)  has  the 
form  expressed  by  Oi .  In  Section  3.3.1,  we  have  shown  that  the  A-S  scheme  can 
satisfactorily  predict  the  net  heating  and  drying  rates  due  to  cumulus  clouds 
for  the  cloud  band  event.  Now,  the  individual  contributions  due  to  the 
condensation,  evaporation , and  vertical  transports  of  heat  and  moisture  to  the 
net  heating  and  drying  rates  are  examined  and  compared  with  Tao's  (1983) 
results  derived  from  the  cumulus  ensemble  model. 


Fig.  3.23  shows  the  distributions  of  heating  effects  due  to  each  component 
at  09  GMT  from  the  semi-prognostic  results.  For  comparison  purposes,  the 
results  were  transformed  into  p-coordinates  and  in  the  units  of  C  hr~^. 


Apparently,  condensation  plays  the  dominant  role  in  the  net  heating  effect. 
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Fig.  3.23.  The  vertical  profiles  of  heating  rate  by  condensation  (C), 
evaporation  (e),  and  vertical  convergence  of  heat  flux  (F).  The  net  heating 
rate  predicted  by  the  A-S  scheme  is  indicated  by  "p",  and  the  net  heating  rate 
etimated  from  large-scale  observations  is  indicated  by  "0". 


Evaporation  and  the  convergence  of  eddy  heat  transport  are  smaller  by  one  order 
of  magnitude.  There  Is  no  evaporation  at  the  low  and  middle  levels  because  no 
cloud  tops  lie  around  those  levels.  Even  though  some  extremely  shallow  clouds 
are  present  (see  Fig.  3.11),  there  is  still  no  evaporation  because  no  detrained 
liquid  water  is  produced  by  the  model  at  that  level.  The  corresponding  results 
simulated  by  Tao  using  the  three-dimensional  version  of  the  cumulus  ensemble 
model  have  been  shown  in  Fig.  3.21.  It  is  seen  that  the  A-S  scheme 
underestimates  the  rate  of  evaporation  by  one  order  of  magnitude.  This  is  not 
surprising  at  all  since  the  A-S  scheme  does  not  include  the  evaporation  of 
rainwater  and  the  evaporation  of  cloud  water  is  assumed  to  take  place  only  at 
the  cloud  tops.  However,  very  interestingly,  the  A-S  scheme  also 
underestimates  the  rate  of  condensation  by  a  factor  of  about  1.6.  Hence,  the 
net  heating  rate  is  about  the  same  as  that  computed  by  Tao  (1983).  Finally, 
the  magnitude  of  the  convergence  of  vertical  eddy  transport  is  generally  much 
smaller  than  the  condensation  in  both  cases. 

To  alleviate  this  deficiency,  the  inclusion  of  convective-scale  downdrafts 
and  their  associated  evaporation  of  falling  precipitation  into  the  scheme  may 
be  desirable.  In  the  next  section  the  results  of  an  experiment  incorporating 
the  downdraft  effects  into  the  A-S  scheme  are  presented. 

3.4  Downdraft  effects; 

We  do  not  intend  to  exaggerate  the  importance  of  downdraft  effects  by 
simply  pointing  out  the  deficiency  encountered  in  the  semi-prognostic  study 
mentioned  above.  Actually,  as  far  as  a  large-scale  prediction  model  is 
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concerned,  the  Importance  of  convective-scale  downdrafts  may  lie  In  the 
interaction  between  the  cumulus  clouds  and  the  boundary  layer  through  the 
detralnment  of  the  downdrafts  Into  the  boundary  layer  and  their  associated 
evaporation  of  the  rainfall  within  the  boundary  layer.  Simpson  and  Van 
Helvolrt  (1980)  discussed  the  Interaction  between  the  clouds  and  sub-cloud 
layer  In  GATE  and  postulated  the  cumulus  downdrafts  to  be  the  main  feedback 
mechanism.  To  translate  this  complicated  physical  process  Into  a 
parameterization  scheme  Is  not  a  straightforward  problem.  First,  we  need  to 
Insure  that  the  cumulus  parameterization  scheme  with  the  Inclusion  of 
downdrafts  will  behave  well.  Secondly,  we  need  to  have  a  boundary  layer 
parameterization  which  can  resolve  the  cumulus  downdraft  effects  well. 

In  this  section,  as  a  first  attempt  to  study  the  downdraft  effects  In  the 
cumulus  parameterization,  the  semi-prognostic  study  is  used  to  examine  whether 
the  accuracy  obtained  in  the  case  without  the  downdrafts  can  be  still 
maintained  and  the  deficiency  mentioned  in  the  preceding  section  can  be 
deleted.  A  mixed  layer  approach  for  the  boundary  layer  parameterization  with 
the  downdraft  effects  will  be  reported  In  the  next  chapter,  where  a  fully 
prognostic  model  Is  developed.  The  discussion  of  the  interaction  between  the 
cumulus  downdrafts  and  the  sub-cloud  layer  will  be  presented  in  Chapter  5. 

3.4.1  Formulation: 

The  downdraft  model  in  this  study  is  adopted  from  Johnson  (1976). 
Basically,  Johnson  proposed  that  each  updraft  cumulus  in  the  original  A-S 
scheme  be  accompanied  by  a  downdraft  element,  which  is  modeled  as  an  "inverted" 
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one-dimens  lonal ,  saturated,  entraining  plume  that  entrains  the  same 
environmental  air  as  the  updraft  does.  The  physical  interaction  between  the 
two  drafts  is  neglected,  except  it  is  assumed  that  the  updraft  supplies  the 
precipitation  to  the  downdraft.  Each  downdraft  is  assumed  to  originate  at  a 
level  given  by  three-fourths  of  the  pressure-depth  of  the  updraft  and  to 
detrain  only  below  the  cloud  base.  The  model  is  closed  by  specifying  the 
relative  intensities  between  the  updraft  and  the  downdraft.  Fig.  3.24 
Illustrates  the  fundamental  features  of  this  model.  Generally  speaking,  the 
nature  of  this  model  is  rather  crude.  However,  in  view  of  that  the  original 
A-S  scheme  is  already  so  complicated,  the  crudity  of  the  downdraft  model  may 
have  its  current  necessity. 


The  details  of  applying  this  downdraft  model  to  the  A-S  scheme  for 
diagnostic  purposes  have  been  reported  in  Johnson  (1976,  1980).  Here,  we  shall 
consider  its  application  to  prognostic  studies.  The  counterparts  of  Eqs. 
(2.1)  and  (2.2)  for  the  downdraft  case  are  now  expressed  by: 
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where  subscripts  u  and  d  denote  updraft  and  downdraft,  respectively,  e^^  is  the 
evaporation  rate  of  liquid  water  detrained  from  the  updrafts,  is  the 
evaporation  rate  of  falling  precipitation  in  the  downdrafts,  and  C  is  the 
condensation  rate,  which  only  occurs  in  updrafts.  Using  the  one-dimensional 
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Fig.  3.24.  The  schematic  diagram  for  the  updraft  and  downdraft  of  cloud  type 


entraining  cloud  model  for  both  the  updrafts  and  downdrafts,  the  vertical  eddy 
fluxes  pw's'  and  pw'q'  now  include  two  parts  expressed  by  (cf.  Eqs.  (2.3)  and 
(2.4)): 


p  w'o'  =  {m^(z  ,z)(a^^j(z,z)  -  a(z) 


+  mj(z,z)(a^jj(z,z)  -  o(z)}dz 


(3.12) 


where  o  =  s,  q,  and  h,  and  Op^j  are  the  thermodynamic  properties  of  the 
updraft  and  downdraft,  respectively.  The  procedure  of  determining  and  Op^j 
is  presented  in  Appendix  A.  The  counterparts  of  (2.5)  and  (2.6)  are 
respectively  written  as: 
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were  2.q{z'>  is  the  downdraft  originating  level  (see  Fig.  3.24).  the  net 

cloud  mass  flux,  is  given  by 


Me  =  Md  + 


(3.18) 


As  mentioned  above,  in  order  to  close  the  model,  Johnson  proposed  a  function 

A 

e(z)  specified  by: 


mQ(z) 
mo(z ) 


(3.19) 


From  observational  indications,  Johnson  proposed  e  =  -0.2  to  -0.4.  It  should 
be  noticed  that  the  parameter  e  does  not  imply  any  significance  of  the 
interaction  between  the  updraft  and  downdraft  but  merely  serves  as  a  closure 
parameter.  Respectively  substituting  (3.13a)  and  (3.14a)  into  (3.10)  and 
(3.11)  with  e^  =  -D4  as  before,  we  obtain  the  counterparts  of  (2.7)  and  (2.8): 


p  +  p^.Vi  +  PW  = 


-  DU  +  (M„  +  Mh) 


(3.20) 


p  +  pV*7q  +  pw  =  D(q*  -  q  +  £)  +  (M  +  M^j) 

9  C  Q  Z  0  Z 


(3.21) 


Comparing  (3.20)  and  (3.21)  with  (2.7)  and  (2.8),  it  is  observed  that  the  net 

effects  on  the  temporal  changes  of  the  thermodynamic  fields  due  to  the 

downdrafts  are  through  the  cumulus-induced  lifting  processes  Mj  —  and  Mj  — . 

^  “  3z 

Mathematically,  with  the  use  of  (3.19),  (3.20)  and  (3.21)  are  still  two 
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Volterra  Integral  equations  with  two  unknowns  mg(z)  and  £(z)  and  they  can  be 
used  in  diagnostic  studies  (Johnson,  1976,  1978,  and  1980). 


s* 

o 


For  prognostic  purposes,  the  next  step  is  to  apply  the  quasi-equilibrium 
assumption  to  each  cloud  sub-ensemble.  To  do  so,  we  take  the  cloud  work 
function  for  updraft  (A^^(z))  into  consideration  by  defining: 


AuCz)  *  /z  - 2 - nu(z,z)(sj.  (z,z) 

CpT(z) 


-  i(z) ]dz 


(3.22) 


Then,  the  counterparts  of  (2.26)  through  (2.30)  can  be  respectively  written  as: 
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Comparing  (3.23)  to  (3.27)  with  (2.26)  to  (2.30),  the  main  difference  between 

the  two  systems  is  the  second  integral  term  in  (3.26),  which  increases  the 

cloud  work  function  of  the  updraft  with  time  due  to  the  induced  lifting  by 

downdrafts.  It  is  also  noticed  from  the  last  two  terms  in  (3.23)  that,  since 

the  inclusion  of  downdrafts  can  modify  the  thermodynamic  properties  of  the 

3h{^ 

mixed  layer  due  to  the  detrainment  of  downdrafts,  -  has  been  divided  into  two 

3hj^  3hj,, 

parts:  and  (-j-j-)cd*  represents  the  time  change  of  hj^  due  to 

large-scale  horizontal  advection  and  the  turbulent  eddy  fluxes  of  heat  and 

moisture  at  the  surface  and  the  top  of  the  mixed  layer  (e.g.,  rhe  right  hand 

3 

side  of  (2.21)).  Meanwhile,  (— — )  .  represents  the  time  change  of  hw  due  to 

the  detrainment  of  cumulus  downdrafts.  However,  for  the  time  being,  the 
3hj^ 

assumption  that  =•  0,  as  described  in  Section  3.2,  is  still  employed  in  this 

semi-prognostic  approach  as  a  simplification,  thus  making  the  last  two  terms  in 

3  h^^ 

(3.23)  vanish.  We  shall  derive  a  formulation  for  ( — in  Chapter  4,  where  a 

3 1 

mixed  layer  parameterization  is  developed  for  a  fully  prognostic  model. 


Finally,  applying  the  argument  of  the  quasi-equilibrium  assumption  to 


(3.23)  and  using  the  algorithm  described  in  Section  3.2,  we  can  solve  for  m^. 
Then,  using  (3.13a)  and  (3.13b),  we  can  determine  the  condensation  rate  (C)  in 
updrafts  and  the  evaporation  rate  due  to  precipitation  (e^)  in  downdrafts. 

3.4.2  The  results: 


We 

have 

extensively  applied 

the  method 

described 

in 

the  previous  section 

to  the 

cloud 

band  event . 

The 

value  chosen 

for 

the 

ratio  between 

the 

intensities 

of  updraft 

and 

downdraft 

(e) 

is 

-0.3 

.  This  value 

was 

experimentally  determined  based  on  the  comparison  between  the  result  derived 
from  the  semi- prognostic  study  at  09  GMT  and  the  averaged  result  from  cloud 
model  simulation  for  the  period  from  06  to  12  GMT  by  Soong  and  Tao  (1980, 
1984).  Fig.  3.25  shows  the  net  cloud  mass  flux  (M^.;  (cf,  (3.18))  as  a 
function  of  time  and  height.  Comparing  Fig.  3.25  with  Fig.  3.9,  it  is  found 
that  the  general  structure  and  the  magnitudes  of  M^  are  similar  to  those  of  M, 
which  again  indicates  that  the  mass  flux  in  the  cloud-free  area  is  negligibly 
small.  Comparing  Fig.  3.25  with  Fig.  3.10,  which  is  the  M^  distribution  for 
the  no-downdraft  case,  it  is  observed  that  the  inclusion  of  downdrafts  would 
not  largely  Influence  the  magnitudes  of  M^.  Fig.  3.26  shows  the  cloud  base 
mass  flux  (mg)  as  a  function  of  cloud  top  height  and  time.  Comparing  Fig.  3.26 
with  Fig.  3.11,  it  is  found  that  the  general  structures  with  and  without 
downdrafts  are  quite  similar.  However,  the  maximum  of  m^  at  the  lower  levels 
in  Fig.  3.26  is  about  two  times  larger  than  that  in  Fig.  3.11,  while  the 
magnitudes  at  the  higher  levels  of  those  two  diagrams  are  about  the  same.  This 
is  not  surprising,  since  the  highest  downdraf t-originating-level  in  this  model 


2 

(GMT) 


is  about  8  km  and,  thus,  all  the  downdraft  effects  are  confined  below  that 
level.  As  discussed  earlier,  the  downdrafts  can  introduce  a  cooling  effect  to 
the  environment  so  that  the  population  of  shallow  clouds  increases.  Note  that 
the  magnitudes  of  with  and  without  downdrafts  are  quite  close,  as  shown  in 
Figs.  3.25  and  3.10,  simply  because  also  includes  the  component  of 

(cf.  (3.18)  and  (3.16)). 

Figs.  3.27  and  3.28,  respectively,  show  the  profiles  of  predicted  Qj^  -  Qj^ 
and  Q2  in  contrast  to  the  observed  ones.  It  is  found  that  the  agreement 
between  the  predicted  and  observed  -  Qj^  is  generally  as  good  as  that  shown 
in  Fig.  3.7  for  the  no-downdraft  case.  However,  in  Fig.  3.28,  less  agreement 
between  the  predicted  and  observed  Q2  at  the  low  levels  is  found  for  most  cases 
which  underestimate  the  cloud  drying  effects.  These  discrepancies  can  be 

expected  because  the  inclusion  of  downdrafts  increases  shallow  clouds,  as 
described  above,  which  would  introduce  an  extra  moistening  effect  by  the 
entrainment  of  water  vapor  into  the  environment  at  the  cloud  tops  (cf.  Eq. 
(2.8)  and  Fig.  3.13)).  However, it  should  be  noted  that  the  sensitivity  of  this 

moistening  effect  is  also  controlled  by  the  relative  humidity.  Thus,  if  the 

lower  atmosphere  is  extremely  moist,  this  discrepancy  would  not  take  place. 

Fig.  3.29  shows  the  results  of  the  predicted  condensation,  evaporation, 
and  vertical  convergence  of  eddy  heat  flux  in  p-coordinates  at  09  GMT.  With 
the  inclusion  of  downdrafts,  it  is  clearly  seen  that  both  the  condensation  and 
evaporation  are  much  larger  than  those  without  downdrafts,  shown  in  Fig.  3.23, 
and  quite  close  to  the  results  derived  from  the  cloud  model,  shown  in 
Fig.  3.21.  Fig.  3.30  shows  the  mass  fluxes  due  to  the  updrafts  and  downdrafts. 


HEIGHT  (km)  HEIGHT  (km) 


HEIGHT  (km) 


Comparing  Fig.  3.30  with  Fig.  3.22.  which  shows  the  cloud  model  results,  both 

the  updraft  and  downdraft  mass  fluxes  predicted  by  the  A-S  scheme  are  found  to 

be  smaller  roughly  by  25%  for  the  former  and  by  35%  for  the  later  than  those 

shown  in  Fig  3.22.  One  possible  reason  for  this  discrepancy  may  be  due  to  the 

different  definitions  for  the  cloud  area.  The  cloud  area  in  the  cloud  model  is 

—  9  —  1 

defined  as  any  area  having  a  liquid  water  content  of  10  g  Kg  or  more.  This 
definition  Includes  active  and  inactive  clouds,  the  later  including  some 
shallow  stratus  clouds  and  anvil  clouds.  On  the  other  hand,  the  cloud  model 
used  in  the  A-S  scheme  only  considers  active  cumulus  clouds.  Nevertheless,  the 
general  agreement  in  the  shape  of  the  mass  flux  profiles  between  the  two 
results  is  quite  well. 

In  passing,  it  should  be  clarified  that  this  downdraft  model  in  the 
prognostic  applications  is  still  in  an  experimental  stage  simply  because  there 
has  been  no  theoretical  way  to  determine  e.  The  same  method  has  been  applied 
to  the  tropical  composite  waves.  It  was  found  that,  in  order  to  obtain 
reasonable  cloud  heating  and  drying  effects,  the  lower  limits  of  e  used  for 
eight  categories  vary  between  0  and  -0.4.  Thus,  concerning  the  incorporation 
of  this  downdraft  model  into  a  large-scale  prediction  model,  the  problems 
become  two-fold:  what  is  the  lower  limit  of  e  for  a  specific  weather  condition? 
and  what  is  the  optimum  e  between  the  lower  limit  and  0?  First  of  ail,  it  is 
felt  that  similar  studies,  as  described  in  this  section,  with  different  data 
over  different  regions  are  required  before  attempting  to  analytically  determine 


the  optimum  e 


CHAPTER  IV 


THE  MESOSCALE  MODEL 


In  Chapter  3,  we  tested  the  A-S  scheme  using  the  semi-prognostic  approach 
to  examine  the  response  of  cumulus  convection  to  the  destabilization  of  the 
atmosphere  by  large-scale  processes.  The  results  indicate  that  the  predicted 
cloud  mass  flux  and  its  associated  heating  and  drying  effects  by  the  A-S  scheme 
are  quite  realistic,  also  bearing  considerable  resemblance  to  those  obtained 
from  a  cloud  modeling  study.  However,  this  close  agreement  between  the 
predictions  and  observations  in  the  semi-prognostic  approach  does  not 
necessarily  imply  .  that  the  A-S  scheme  will  have  prognostic  value  in  a  fully 
prognostic  model.  This  point  should  be  clear,  since,  in  a  fully  prognostic 
model,  both  the  cloud  properties  and  large-scale  disturbances  evolve  in  time 
and  highly  Interact  with  each  other. 

In  order  to  examine  how  the  A-S  scheme  works  in  a  complete  prognostic 
problem,  a  two-dimensional  hydrostatic  model  was  developed  with  the  inclusion 
of  the  A-S  scheme  to  investigate  the  evolution  and  structure  of  the  same 
tropical  cloud  band,  discussed  previously,  under  the  influence  of  a  low-level 
large-scale  lifting.  This  chapter  will  describe  the  model.  The  model  results 
will  be  represented  in  Chapter  5. 
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4.1  Governing  equations 


Although  mesoscale  dynamics  Is  characteristically  three-dimensional,  as  a 
first  attempt  to  understand  the  response  of  a  tropical  mesoscale  convective 
system  to  the  large-scale  forcing  through  the  application  of  the  A-S  scheme,  we 
assume  that  the  predicted  quantities  do  not  depend  on  one  of  the  horizontal 
Cartesian  coordinates  (namely  y).  In  order  to  Incorporate  rotational  effects, 
we  retain  the  velocity  component  In  the  y  direction.  The  mesoscale  model 
presented  here  employs  the  "deep  anelastlc"  equations  (Ogura  and  Phillips, 
1962)  but  with  the  hydrostatic  assumption.  The  two-dimensional  "hydrostatic 
version"  of  the  anelastlc  system  has  been  used  by  Orlanski  and  Ross  (1977)  with 
the  replacements  of  the  momentum  equations  by  the  vorticity  equation  in  a 
frontal  simulation.  Rosenthal  (1969,  1970)  used  this  version  but  In  an 
axisymmetrlc  model  to  simulate  the  development  of  hurricanes. 


The  momentum  equations  In  the  x-z  slab  symmetric  domain  are: 

(4.1) 

(4.2) 

where  0  is  the  potential  temperature,  f  Is  the  Coriolis  parameter  with  the 
value  of  2.15  x  10”^  s“^,  Ug  is  the  x-component  of  the  geostrophic  wind,  which 
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is  taken  Co  be  Che  constant  value  of  Che  initial  wind  u  throughout  the 
Integration,  and  n ,  the  nondlmensional  pressure,  Is  defined  as; 


(4.3) 


where  p  Is  pressure,  Pq  Is  the  reference  pressure  (1000  mb),  R  is  the  gas 
constant  of  dry  air,  and  is  the  specific  heat  of  dry  air  at  constant 
pressure.  The  hydrostatic  equation  is: 


'(■“v  57  ■  -8 


(4.4) 


where  g  is  the  gravity  and  0  is  the  virtual  potential  temperature  defined  by; 


0y  =>  0(1  +  O.eiq) 


(4.5) 


The  thermodynamic  and  moisture  equations  are: 
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where  is  the  dlabatlc  heating  rate  due  to  clouds,  is  the  production  rate 

96 

of  water  vapor  due  to  clouds,  and  _ 8  related  to  u„  by  the  thermal  wind 

3y  ^ 

relation: 
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3y  g  32 


(A.8) 


where  6„  can  be  specified  from  the  initial  soundings.  It  is  noted  that 

O 

radiation  effects  are  neglected  in  (4.6),  assuming  that  they  are  not 
fundamental  to  the  development  of  a  tropical  mesoscale  convective  system.  The 
mass  continuity  equation  is  expressed  in  the  anelastic  form: 


iii  +  i  =  0 

9x  p  dz 


(4.9) 


where  p,  the  density,  is  a  function  of  height  only. 


Respectively  separating  the  variables  0,  q,  and  it  into  the  horizontal  mean 
and  its  departure  as : 


A(x,z,t)  =  A(z,t)  +  A'(x,z,t) 


(4.10) 


where  A  =*  9  ,  q,  and  it,  and  subsequently  substituting  (4.10)  into  (4.1),  (4.2), 
(4.4),  (4.6),  and  (4.7),  with  the  use  of  (4.9),  we  obtain  a  set  of  governing 
equatl'^ns  in  flux  form: 


3u  3uu  13,  \.c  Z 

—  - - - (puw)  +  fv  -  c„9  - 

3t  3x  p  3z  P  3x 


(4.11) 


3v  3uv  13,  X  d  \ 

3 1  3x  p  3z  8 


(4.12) 
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Cp0  ^  -  g(il  +  0.61q') 
0 


(A. 13) 


3t 


^  (u0')  - 
3x  p  3z 


(pwe  ' )  -  V  -  w 

3y  3z  '- 


(A. 14) 


?Ji  =  -  i_ 

3t  3x 


(uq') 


i  _L 

p  3z 


(pwq')  -  w  +  P 
dZ  ' 


(A. 15) 


In  deriving  (4.13),  we  assumed  chat  the  horizontal  means  satisfy  the 
hydrostatic  relation; 


Vvfl'-g  (A.  16) 

In  deriving  (4.11),  0  in  front  of  the  pressure  gradient  term  was  approximated 
by  0 . 

Since  this  model  is  designed  to  examine  the  evolution  of  a  mesoscale 
convective  system  under  the  Influence  of  a  weak  large-scale  lifting,  a 
constantly  imposed  low-level  circulation  [wq,  uq) ,  which  satisfies  the 
continuity  equation,  is  now  introduced  into  (4.11),  (4.12),  (4.14),  and  (4.15): 

It  *  ■  (u+uo)(u+Uo)  -  (p{u+u^j)(w+Wq)] 

+  fv  -  C.0  ^ 

P  3x 


(A. 17) 


39  3  /  ^  \a'  13,^,.,  ®“g 

_.--^„«0)9  -__(«vo)6  -.— 

-<“«o)  5T  +  Ol  <*■'’> 


If  ■  ■  H  '“■“o)"'  - 1  ^  ('«'o)i'  -  (-"^o)  If  +  fv 


(4.20) 


The  specifications  for  Wg  and  Ug  are  discussed  In  Section  4,3,2. 


The  final  step  is  to  replace  the  cloud  terms  represented  by  and  in 
(4,19)  and  (4.20)  by  the  A-S  parameterization.  Using  the  notations  introduced 
in  Chapter  2,  we  have: 


=  _L_  [Me  II  -  DU] 


(4.21) 


Since  the  thermodynamic  equation  in  the  A-S  scheme  was  originally  formulated  in 
terms  of  the  static  energy  (s),  we  need  to  convert  0  into  s,  thus  requiring  thq 
use  of  the  nondimensional  pressure  n.  In  doing  this,  we  adopt  the 


approximation  that  tt'/w  is  an  order  of  magnitude  less  than  9'/^,  from 
Wllhelmson  and  Ogura  (1972): 

—  «  —  (4.23) 

iT  9 

such  that  T  Is  approximated  by  T  *  ^9  and  s  is  approximated  by  s  *  CpiF9  +  gz. 

Now,  the  set  of  equations  which  govern  the  evolution  of  a  mesoscale 
convective  system  under  the  Influence  of  the  large-scale  lifting  consists  of 
(4.17)  to  (4.20)  with  the  continuity  equation  (4.9)  and  hydrostatic  equation 
(4.13). 

4.2  Boundary  layer  parameterization: 

A  mixed  layer  approach  (see  Fig.  2.1)  is  employed  with  the  assumption  that 
the  depth  of  the  mixed  layer  is  constant  in  time  and  space.  The  equations  for 
and  q]^^  are  written  as; 
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(4.24) 
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3z 
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(4.25) 


where  pw'9 '  and  pw'q'  represent  the  cumulus  fluxes  of  heat  and  moisture,  while 


where 


Az  “  z^,  subscript  B  denotes  the 
surface.  It  Is  assumed  that  the  eddy 
and  the  cumulus  fluxes  at  the  surface 

(^u>B  =  (Vb  =  0 


cloud  base,  and  subscript  0  denotes  the 
momentum  fluxes  at  the  top  of  mixed  layer 
vanish: 

(4. 34) 


(ow'0  ')q  »  (pw'q'  )g 


(4.35) 


Using  Deardoff's  parameterization  (1972),  we  can  determine  (f0)o»  (^q^O* 
(Fv)0  by: 


(F«) 


9^0 


9v  -®v 
'^0 


(4.36) 


(qo“‘iM)(Fe  ) 
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(4.38) 


where  (Fg)Q  is  the  eddy  flux  of  virtual  potential  temperature  at  the  surface 
and  Is  given  by: 


v' 

•  ' 

E 


a 


-  »0<=9“.«.o  -  «.„> 


(4.39) 


The 

notations 

Introduced 

in  (4.36) 

to 

(4.39)  are  explained 

as  follows : 

Pq 

the 

air  density  at  the 

surface. 

‘^u 

the  friction  coefficient,  Cg 

the 

heat¬ 

-transfer 

coefficient, 

and  u^ 

the 

frictional  velocity. 

u«  Is  determined 

hy: 

i2 

PQ 


(4.40) 


and  Cg  are  determined  through  empirical  formulas  which  are  functions  of  the 
bulk  Richardson  number 


,  .  s 
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(4.41) 


Following  Lilly  (1968),  we  determine  the  turbulent  eddy  fluxes  of  heat  and 
moisture  at  the  top  of  the  mixed  layer  by: 


(Fg)B  -  -  k(F„) 


e-'o 


(4.42) 


(4.43) 


•V  ‘,  *  T 


mm.*  ^  .V  .  \ 


::4 


where  k  Is  chosen  to  be  0.2. 


The  next  step  Is  to  determine  the  cumulus  fluxes  of  heat  and  moisture  at 

the  top  of  the  mixed  layer,  which  are  respectively  represented  by  (pw'0')g  and 

(pw'q')g  in  (4.28)  and  (4.29).  If  following  the  Arakawa  and  Schubert's 

approach  as  stated  in  Sections  2.1,  there  would  be  no  contributions  to  (pw'0')g 
and  (pw'qjg  due  to  the  updrafts.  This  can  be  shown  by  the  relations  (cf. 

(2.3))  and  (2.4)); 

(pw^0')g  -  /2'“*^nB(2')nu(zg,z')(0j.^(zg,z')  -  0(zg))di  (4.44) 

u  B 

(pw'q')B  =  ~  q(zg))dz  (4.45) 

u  B 

A 

where  subscript  u  represents  the  updraft,  and  with  9cu(zb»z)  “ 
q^,y(zB.2)  =  q(zB^  obtain  (pw'0  ')gy  =  (pw'q)gy  =  0.  Hence,  In 

this  model,  (ow'0')g  and  (pw'q')g  are  only  associated  with  downdrafts  and  can 
be  formulated  by: 

(pw'o')g  =  j  2^^  ~  o(zB+)))dz'  (4.46) 

where  a  -  0  and  q.  mQ(z)  Is  the  downdraft  mass  flux  of  sub-ensemble  z  at  Its 
downdraft  originating  level  and  njj(zg,z)  is  the  normalized  mass  flux  for 
downdraft.  thermodynamic  property  of  the  cumulus  downdraft  of 

sub-ensemble  z  at  the  top  of  the  mixed  layer  and  a(zg+)  is  the  thermodynamic 
property  of  the  environment  at  the  level  "immediately"  above  the  top  of  the  the 
mixed  layer,  assuming  that  an  infinitesimally  thin  transition  layer  is  above 
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the  top  of  the  the  mixed  layer  (see  Fig.  2.1).  Referring  to  Fig.  3.4,  the 
value  of  a  at  2^+  Is  the  one  interpolated  from  levels  I  and  2.  On  the  other 
hand ,  the  value  of  a  at  Zg  is  the  same  as  the  one  for  level  1 .  This  treatment 
of  the  cumulus  downdraft  fluxes  at  the  top  of  the  mixed  layer  is  essentially 
the  same  as  Johnson's  (1976). 

Applying  (4.46)  to  the  moist  static  energy  (h),  we  have: 

(p^^)b  -  /f®®*  mQ(z')n^i(zB,^)(h^^(2B,^)  -  h(2B+))d^  (4.47) 

B 


Now,  referriag  back  to  (3«23),  we  can  express  ( 


ay 

TT^cd 


by : 


a  ^ 

^IT^cd 


1  /*max 


(4.48) 


-  h(2g+))dz' 

With  the  uses  of  (4.48)  and  (3.19),  the  last  term  In  (3.23)  can  be  expressed 
as : 


3hM  ! 

<lf>cd  !; 


(4.49) 


.fz  p(z)0(z)dz  (4.50) 

Namely,  we  have  Introduced  an  extra  kernel  into  (3.23),  which  arises  due  to 
the  presence  of  Che  downdrafts. 

To  determine  the  evaporation  rate  of  falling  precipitation  in  the 
sub-cloud  layer  (e^^  (4.28)  and  (4.29),  we  simply  take  the  Interpolated 

value  of  e^  (see  Eqs.  (3.13c)  and  (3.14c))  from  level  1  (see  Fig.  A. 3)  and  the 
surface,  where  we  assume  that  e^  vanishes. 

4.3  Model  design: 

4.3.1  Model  structure  and  finite  difference  schemes: 

As  shown  in  Fig.  4.1,  this  model  has  a  horizontal  domain  of  1600  km,  with 
600  km  of  fine  grids  at  the  center  part  of  the  domain  with  a  grid  spacing  of 
20  km.  On  each  side  of  the  fine-grid  area,  there  is  an  environmental  area  with 
10  grid  points  stretched  by  a  ratio  of  1.2  to  the  preceding  grid  size.  17  grid 
points  are  located  in  the  vertical  z-coordinates  with  an  even  grid  spacing  of 
900  m  such  chat  the  top  of  the  domain  is  at  the  level  of  14.4  km,  which  is 
around  140  mb  in  p-coordinates .  The  mixed  layer  parameterization  described  in 
Section  4.2  is  applied  to  the  first  layer  of  the  model  such  that  the  cloud  base 
height  is  fixed  at  900  m.  Since  the  convective  system  under  investigation  is 


-I  co¬ 


in  the  form  of  a  line  structure,  as  shown  in  Fig.  3.5,  the  x-coordinates  in  the 
model  are  chosen  in  the  direction  normal  to  the  cloud  band  with  the  positive  x 
direction  toward  the  NNW  (Soong  and  Tao,  1980). 

A  leapfrog  time  integration  and  a  second-order  centered  space  difference 
are  applied  with  a  staggered  grid  arrangement,  shown  in  Fig.  4.2.  The  finite 
difference  equations  corresponding  to  (4.9)  and  (4.17)  through  (4.20)  are  given 
in  Appendix  C.  A  time-smoothing  scheme  formulated  by  Robert  (1966)  is  adopted 
in  the  integration  to  avoid  the  time  splitting.  The  scheme  is  of  the  form: 

u*"^^  =  u''"^  +  F"  2At  (4.51a) 


u”  =  u 


*n 


+  0(u 


-  2u 


*n 


u"“^) 


(4.51b) 


where  superscript  n  denotes  the  time  step,  *  denotes  the  unsmoothed  value,  F 
represents  all  the  terms  on  the  right  hand  side  of  the  prognostic  equation,  and 
a  is  the  smooching  coefficient.  A  value  of  0.1  is  used  for  a.  The  model  also 
includes  a  horizontal  fourth-order  numerical  smoother  (Klemp  and  Wilhelmson, 
1978;  Bradley,  1985)  which  is  applied  to  all  of  the  prognostic  fields.  The 
smoother  is  of  the  form: 

“  u(^i+2  '^i-2  ”  ^^'^1+1  '^i-P  ^*^i^  (^*52) 

where  the  coefficient  u  is  0.035.  This  horizontal  smoother  reduces  the  growth 
of  nonlinear  instabilities  and  filters  out  very  short  waves.  The  use  of  this 
fourth-order  smoother  makes  Che  damping  effects  very  selective  by  attenuating 
the  short  waves  very  effectively  but  leaving  the  longer  waves  almost  unchanged. 


I 

We  apply  a  second-order  smoother  to  Chose  grid  points  adjacent  to  Che  lateral 
I  boundaries.  In  this  model,  no  vertical  smoothing  Is  Included.  In  addition  to 

satisfying  the  Courant-Frledrlchs-Levy  (CFL)  condition,  the  determination  of 
the  time  step  (At)  in  this  kind  of  primitive  mesoscale  model  requires  caution, 
because  of  the  behavior  of  gravity  waves  associated  with  the  hydrostatic 
assumption  (Tapp  and  White,  1976).  In  a  non-hydrostatic  atmosphere,  the 
frequency  of  gravity  waves  is  bounded  by  the  Brunt-Vaisala  frequency  N. 
However,  in  a  hydrostatic  atmosphere,  the  frequency  of  the  gravity  waves  goes 
to  infinity  as  the  wave  number  increases.  Thus,  all  of  the  gravity  waves  with 
frequencies  larger  than  N  in  a  hydrostatic  atmosphere  should  be  regarded  as 
high-frequency  spurious  waves  that  can  introduce  numerical  instability  if  the 
choice  of  AC  is  too  large.  Through  the  "local  instability"  analysis,  the  CFL 
condition  for  this  model  is: 


At  C 


1 


+  I 


Az 


(4.53) 


From  (4.53),  the  magnitude  of  At  can  be  as  large  as  10  seconds.  However,  for 
the  reason  stated  above,  we  choose  At  to  be  90  seconds  for  the  model 
integration. 

4.3.2  Boundary  and  initial  conditions: 

Fixed  upper  and  lower  boundary  conditions  are  applied  for  w,  while  zero- 


gradient  lateral  boundary  conditions  are  applied  for  9  and  q  and  open  lateral 
boundary  conditions  for  u  and  v: 
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w  »  0 


at  z  ■  0,  z>^ 


(4.54) 
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(4.55) 


where  z.j.  represents  the  top  of  the  domain,  -£  and  £  represent  the  left  and 
right  boundaries,  respectively,  and  c  represents  a  gravity-wave  speed.  A 

O 

value  of  10  m  s“^  for  Cg  Is  used. 

The  soundings  (shown  In  Fig.  4.3)  taken  from  the  ship  Researcher  at  00  GMT 
12  August  1974  are  used  for  the  background  Initial  conditions  for  9  and  q. 
Fig.  4.4  shows  the  low-level  vertical  velocity  at  00  GMT  12  August  1974 
computed  from  the  wind  data  over  the  GATE  A/B  array.  It  Is  seen  that  the 
maximum  magnitude  Is  about  1.7  cm  s~^  around  850  mb.  We  double  the  magnitudes 
shown  in  Fig.  4.4  and  multiply  them  by  a  sinusoidal  function  with  the  value  of 
unity  at  the  center  of  the  domain  and  zero  at  the  boundaries  of  the  fine-grid 
area  (see  Fig.  4.1)  to  obtain  the  imposed  vertical  velocity  Wq(x,z)  included  in 
(4.17)  to  (4.20).  The  corresponding  Uq(x,z)  is  calculated  from  the  continuity 
equation  with  the  specification  of  Uq  »  0  at  the  center  line  of  the  model 
domain.  This  [wq,  uqJ  will  be  constantly  imposed  on  the  model  as  a  time- 
independent,  large-scale  forcing.  Fig.  4.5,  adapted  from  Ogura  aT .  (1979), 

shows  the  surface  wind  and  divergence  calculated  at  the  initial  stage  of  the 


PRESSURE,  mb 


Fig.  4.3.  The  atmospheric  soundings  taken  f 
0000  GMT  12  August  1974. 


PRESSURE  (mb) 


lC 


Fig.  A. 4.  The  computed  low-level  vertical  velocity  at  0000  GMT  12 
August  1974  at  the  ship  Researcher. 


LATITUDE 


LONGITUDE 


Fig.  4.5.  The  surface  wind  and  divergence  (in  the  units  of  10~^  s”^) 
along  with  radar  echoes  observed  at  0000  GMT  12  August  1974.  The 
outer  hexagon  indicates  the  GATE  A/B  scale  ship  array  and  the  inner 
hexagon  indicates  B  scale  ship  array.  The  full  barb  for  wind  is 
5  m  s“^  (Adapted  from  Ogura  et  al.,  1979). 
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cloud  band.  It  Is  evident  that  convective  cells  started  developing  In  the  area 
of  the  maximum  surface  convergence.  Almost  symmetric  decreases  of  surface 
convergence  are  present  on  both  sides  of  the  cloud  band.  This  indicates  that 
our  Imposed  circulation  [wq,  uqI  described  above  is  fairly  Justifiable.  The 
data  analysis  also  shows  that  a  700  mb  easterly  wave  was  passing  through  the 
area  during  the  period  considered  and  that  the  A/B  array  was  in  category  3  when 
the  development  of  the  cloud  band  occurred.  This  Implies  that  the  large-scale, 
low-level  lifting  mechanism  Is  associated  with  the  passage  of  the  easterly 
wave. 

4.3.3  Computational  procedure: 

A  computational  cycle  Is  outlined  as  follows: 

(1)  Compute  the  advectlon  and  Coriolis  terms  In  momentum  equations  (4.17) 
and  (4.18). 

(11)  Compute  the  advectlon  terms  In  thermodynamic  and  moisture  equations 
(4.19)  and  (4.20). 

(ill)  Use  the  A-S  scheme  to  compute  the  cloud  effects.  Since  using  the  A-S 
scheme  requires  considerable  computing  time,  the  subroutine  for  the 
A-S  scheme  is  called  every  15  time  steps.  The  values  of  the  cloud 
effects  within  each  13  time  steps  are  the  ones  calculated  at  the 
latest  call  of  the  A-S  scheme.  Once  the  cloud  effects  In  (4.21)  and 
(4.22)  are  determined,  we  can  obtain  6  and  q  for  the  next  time  step. 

(Iv)  Use  the  boundary  condition  (4.54)  and  the  mass  continuity  equation 


(4.9)  to  obtain  the  relation 


Then,  Inserting  (4.17)  Into  (4.56)  and  using  the  hydrostatic  equation 
(4.13),  we  can  diagnostically  determine 

dx 

* 

(v)  Use  - —  obtained  from  (Iv)  and  the  values  obtained  from  (1)  to  advance 
3x 

u  and  V  for  the  next  time  step. 

(vl)  Apply  the  horizontal  smoother  to  8 ,  q,  u,  and  v. 

(vll)  Use  u  calculated  from  (vi)  and  the  continuity  equation  to  determine  w 


for  the  next  time  step. 


CHAPTER  V 


MODEL  RESULTS  OF  MESOSCALE  SIMULATION 


In  this  chapter,  the  results  of  the  response  of  the  tropical  mesoscale 
cloud  band  to  the  large-scale  low-level  lifting  process  through  the 
applications  of  the  A-S  scheme  are  illustrated  in  terms  of  the  following 
"circulation  descriptors"  defined  by  Kreltzberg  and  Perlcey  (1977):  (1)  the 
Intensity  of  the  mesoscale  vertical  motion,  (11)  the  time  scale  of  the 
development  of  the  mesoscale  circulation,  (ill)  Che  time  scale  of  Che  decay  of 
Che  mesoscale  circulation,  (Iv)  the  amount  of  surface  precipitation,  and 
(v)  Che  horizontal  width  scale  of  the  cloud  band  and  Its  associated 
precipitation  area.  For  verification  purposes,  the  model  results  will  be 
compared  with  the  observed  values  deduced  from  data  analysis.  Furthermore,  a 
detailed  examination  of  the  heating  profiles  under  different  mesoscale 
environmental  conditions  (primarily  vertical  velocity  profiles)  will  be  made  to 
investigate  the  roles  played  by  the  convective  heating,  as  represented  by  the 
A-S  scheme.  In  the  evolution  of  the  mesoscale  circulation. 

Based  on  Che  model  characteristics  described  In  Chapter  4,  the  results  of 
5  runs  are  presented  In  this  chapter.  Run  1,  which  Is  hereafter  referred  to  as 
the  control  run,  excludes  the  Initial  horizontal  wind  (u  )  (see  Eq.  (4.2))  and 

O 

the  cumulus  downdraft  effects  in  the  A-S  scheme.  The  results  of  a  24-hour 
Integration  for  the  control  run  are  discussed  in  Section  5.1.  The  other  4 


sensitivity  test  runs  are  presented  in  Section  3.2.  Note  that  only  one  factor 


is  changed  In  each  subsequent  run.  In  run  2,  we  performed  an  experiment  by 
deleting  the  large-scale  low-level,  forcing  [wq,  uq]  at  the  time  when  the 
mesoscale  circulation  reaches  its  peak  Intensity  in  the  control  run.  This 
experiment  would  reveal  whether  or  not  a  well-defined  mesoscale  circulation  can 
sustain  itself  for  a  considerable  period  of  time  without  the  large-scale,  low- 
level  lifting  process.  Run  3  presents  an  experiment  with  the  Inclusion  of  the 
initial  sheared  horizontal  wind  (u^)  since,  based  on  various  observational  and 
simulation  studies,  the  vertical  wind  shear  has  been  considered  very  Important 
to  the  development  and  movement  of  the  mesoscale  convective  system.  Run  4 
presents  an  experiment  with  the  inclusion  of  the  cumulus  downdraft  effects  into 
the  A-S  scheme.  Emphasis  here  will  be  placed  on  the  feedback  mechanism  between 
the  cloud  layer  and  the  mixed  layer  due  to  the  detrainment  of  the  downdrafts. 
Finally,  in  run  5  we  replaced  the  A-S  scheme  by  Kuo's  1974  scheme  for 
comparison  purposes. 

5.1  The  control  simulation: 

5.1.1  Results  of  an  early  simulation: 

The  initial  conditions  shown  in  Fig.  4.3  are  conditionally  unstable  except 
for  a  stable  layer  around  700  mb.  There  is  no  initial  motion  field  except  for 
the  externally-imposed,  time-independent,  large-scale,  low-level  circulation 
[uq,  Wq] ,  as  described  in  Section  4.3.2.  The  Wq  field  is  shown  in  Fig.  5.1. 
Since  the  model  is  symmetric  about  the  center  axis  of  the  domain  (see 
Fig.  4.1),  hereafter  all  of  the  model  results  are  presented  only  in  the  right 
half  of  the  domain  except  for  the  run  with  vertical  wind  shear. 
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The  mesoscale  fields  will  evolve  from  the  above  Initial  conditions  due  to 
the  horlzontally-dlf ferentlal  convective  heating  rates  resulting  from  the 
applications  of  the  A-S  scheme.  Fig.  5. 2a  shows  the  convective  heating  rates 
at  1  hr.  Note  that  since  the  cloud  base  height  Is  set  at  900  m,  no  convective 
heating  Is  observed  below  that  level.  It  Is  clearly  seen  that  significant 
convective  heating  concentrates  In  the  central  portion  of  the  domain  with  a 
horizontally  well-defined  differential  distribution  because  of  the  distribution 
of  Wq,  as  shown  in  Fig.  5.1.  Fig.  5.2b  shows  the  model-predicted  vertical 
velocity  field  (w)  at  1  hr.  Obviously,  the  concentration  of  w  in  the  center 
part  of  the  domain  Is  due  to  the  distribution  of  convective  heating  shown  in 
Fig.  5.2a.  The  maximum  of  w  is  about  2  cm  s~^  and  occurs  at  the  center  and 
around  6  km  in  height.  A  broad  area  of  weak  downward  motion  Is  present  away 
from  the  center.  The  deviation  of  the  potential  temperature  (6')  from  the 
initial  state  is  shown  In  Fig.  5.2c.  The  maxima  of  positive  anomalies  located 
at  5  and  10  km  In  height  at  the  center  columns  are  obviously  due  to  the  strong 
convective  heating.  Negative  anomalies  are  found  at  the  low  levels  around  the 
center  because  of  the  imposed  Wq,  which  results  In  relatively  stronger 
adiabatic  cooling  than  the  convective  heating  derived  from  the  A-S  scheme  for 
the  low  levels.  On  the  other  hand,  the  weak  positive  anomalies  between  120  and 
200  km  are  caused  by  the  adiabatic  warming  due  to  the  downward  motion.  The 
temperature  distributions  shown  In  Fig.  5.2c  are  expected  to  generate  a 
low-level  mesolow  around  the  central  portion  of  the  domain  and  to  cause  a 
circulation  with  the  patterns  of  w  and  u,  as  shown  in  Figs.  5.2b  and  5.3a, 
respectively.  The  fields  of  the  ageostrophlc  wind  (v^g)  defined  by: 


s 


HEIGHT  (Km)  HEIGHT  (Km) 


a 


Fig.  5.3.  Cross  sections 
(3)  x-cooponent  velocity 


(a)  x-component  velocity  In  units  of  cm  s'^  with  contour 
b)  y-component  ageostrophlc  wind  in  units  of  cm  with  contour 
100,  and  (c)  y-component  velocity  in  units  of  cm  s"^  with  conto 


simulation  of 
units  of  cm 


the  control  run  at  1 


hr: 

with  contour  interval  of  7, 
interval  of 
contour  interval  of 


and  the  y-component  velocity  (v)  are  shown  in  Figs.  5.3b  and  5.3c.  Vnierever  v 

Is  sub-geostrophlc  (v__  <  0) ,  the  air  will  accelerate  to  the  left  and  cause 

ag 

convergence  toward  the  center.  Fig.  S.Aa  shows  the  cloud  base  mass  £1  x  (mg) 
distribution  as  a  function  of  cloud  top  height  and  distance  from  the  center  at 
1  hr.  It  is  evident  that  the  blmodal  distributions  with  relatively  weaker  high 
clouds  are  present  around  the  central  columns,  while  only  shallow  clouds  are 

present  in  the  outer  part  of  the  domain.  The  corresponding  cloud  mass  flux 

(M^)  is  shown  in  Fig.  5. 4b.  The  total  mass  flux  (M)  is  now  defined  by: 

M  -  p(w  +  Wq)  (5.2) 

and  shown  in  Fig.  S.4c.  It  is  again  found  that  the  magnitudes  of  and  M  are 
very  close  to  each  other,  as  Indicated  in  the  semi-prognostic  studies. 

The  results  of  the  early  simulation  at  1  hr  have  been  discussed  in  the 
preceding  paragraph.  Since  all  of  the  fields  evolving  from  the  initial 
conditions  are  driven  by  the  horizontally  differential  convective  heating, 
which  is  computed  from  the  A-S  scheme,  it  is  important  to  know  how  the  A-S 
scheme  responds  to  the  imposed  large-scale,  low-level  forcing  at  the  very 
beginning.  Fig.  5.5  shows  the  profiles  of  the  adiabatic  cooling  rate  due  to 
the  Imposed  Wq  and  the  convective  heating  rate  deduced  from  the  A-S  scheme  for 
the  central  column.  It  can  be  seen  that  the  adiabatic  cooling  ceases  at  500  mb 
because  Wq  vanishes  above  that  level  while  the  convective  heating  takes  place 
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Fig.  5.4.  Cross  sections  of  simulation  of  the  control  run  at  I  hr:  (a)  the 
cloud  base  mass  flux  (mg)  in  units  of  mb  hr“^  with  the  contour  interval  of  0.6, 
(b)  the  cloud  mass  flux  (M^,)  in  unit^s  of  mb  hr”^  with  the  contour  Interval  of 
2.0,  and  (c)  the  total  mass  flux  (M)  In  units  of  mb  hr~^  with  contour  interval 
of  2.0. 
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from  the  cloud  base  up  to  200  mb.  The  small  differences  between  the  two 
profiles  from  825  mb  to  200  mb  Introduce  positive  temperature  tendencies  and 
the  small  differences  below  825  mb  Introduce  negative  tendencies.  Meanwhile, 
the  differences  between  the  adiabatic  cooling  and  convective  heating  rates 
gradually  decrease  outward  from  the  center  because  of  the  distribution  of  Wq  so 
that  a  horizontal  temperature  gradient  Is  generated.  Note  that,  although  the 
imposed  low-level  forcing  Is  confined  below  550  mb,  high  clouds  are  still 
present  because  the  large-scale  forcing  for  the  clouds  with  cloud  tops  higher 
than  550  mb  does  exist  (see  2.27)). 


5.1.2  A  24-hour  Integration 


Fig.  5.6  illustrates  the  evolution  of  the  mesoscale  circulation  In  the 
vector  fields  of  (u,  w].  The  magnitudes  of  w  have  been  enlarged  by  10  times 
because  of  the  scale  difference  between  u  and  w,  at  3,  9,  and  17  hr,  which  can 
be  respectively  representative  of  developing,  mature,  and  decaying  stages  based 
on  their  flow  structures  and  Intensity  of  vertical  motion.  The  flow  pattern, 
with  inflows  at  low  levels,  upward  motion  at  the  central  portion  of  the  domain, 
and  outflows  at  upper  levels.  Is  distinctly  present  In  Figs.  5.6a  and  5.6b. 
Fig.  5.7  shows  the  evolution  of  the  vertical  velocity  at  the  different  stages. 
It  is  seen  that  at  3  hr  the  magnitudes  of  w  become  much  larger  Chan  chose  shown 
in  Fig.  5.2b.  However,  the  upward  motion  still  concentrates  around  the  center. 
Significant  downward  motion  with  magnitudes  of  about  2  to  3  cm  s”^  Is  present 
beside  the  upward  motion  area.  At  9  hr,  the  Intensity  of  the  upward  motion 
reaches  a  magnitude  of  17  cm  s~^  and  both  the  upward  motion  area  and  the 
significant  downward  motion  area  are  more  broad.  At  17  hr,  the  system  Is 
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5.6.  Cross  sections  of  vector  field  of  [u,w]  in  the  control  run:  (a)  at  3 
(b)  at  9  hr,  and  (c)  at  17  hr.  The  magnitudes  of  w  have  been  enlarged  by 
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Fig.  5.7.  Cross  sections  of  vertical  velocity  (w)  In  units  of  cm  s”  in  the 
control  run:  (a)  at  3  hr  with  contour  interval  of  2.0,  (b)  at  9  hr  with  contour 
interval  of  2.0,  and  (c)  at  17  hr  with  contour  interval  of  1.0. 


120 


240  360  420  600 

HORIZONTAL  DISTANCE  (Km) 


720 


fig.  5.7.  (Continued) 
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obviously  in  Its  decaying  stage  with  strong  downward  motion  located  at  low 
levels  around  the  center  and  weak  upward  motion  at  high  levels  away  from  the 
center.  Note  that  all  of  the  w  fields  presented  in  this  chapter  do  not  include 
the  imposed  Wq.  The  decaying  situation  is  also  indicated  by  Fig.  5.6c  where  a 
well-defined  reverse  circulation  is  present  at  the  low  levels.  Fig.  5.8  shows 
the  field  of  0 ' .  Examining  the  w  field  at  3  hr  shown  in  Fig.  5.7a,  It  should 
be  certain  that  the  positive  temperature  anomalies  at  3  hr  in  the  area  between 
120  to  400  km  are  caused  by  the  adiabatic  warming  due  to  the  downward  motion. 
It  can  be  seen  in  Fig.  5.8b  that  after  the  9-hour  simulation  the  temperature 
Increases  about  1.9  K  for  the  atmosphere  at  11  km  in  height  around  the  center. 
Meanwhile,  the  temperature  decreases  about  0.2  K  for  the  atmosphere  at  low 
levels  around  the  center.  These  features  can  be  readily  explained  by 
Fig.  5.9a,  which  shows  the  profiles  of  the  adiabatic  cooling  rate  due  to 
vertical  and  horizontal  advection  and  the  convective  heating  rate  derived  from 
the  A-S  scheme  for  the  center  column  at  6  hr.  As  discussed  earlier,  the  small 
difference  between  these  two  profiles  is  responsible  for  the  temperature 
changes  in  the  convective  area.  The  double  maximum  distribution  of  the 
adiabatic  cooling,  shown  in  Fig.  5.9a,  is  due  to  the  distribution  of  total 
vertical  velocity  (w  +  Wq)  shown  in  Fig.  5.9b.  The  lower  maximum  is  apparently 
caused  by  the  imposed  low-level  Wq.  Generally  speaking,  the  profile  shown  in 
Fig.  5.9b  characterizes  the  situations  of  strong  convection  in  the  tropics. 

The  gradual  accumulation  of  the  adiabatic  warming  due  to  the  downward 
motion,  as  noted  above,  results  in  a  secondary  maximum  of  the  deviation  of 
potential  temperature  about  300  km  from  the  center  and  11  km  in  height  at  9  hr, 
as  shown  in  Fig.  5.8b.  Of  great  Importance  is  that  this  secondary  maximum  of 
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Fig.  5.8.  Cross  sections  of  the  deviation  of  potential  temperature  (6')  form 
Che  initial  condition  in  units  of  K  in  the  control  run:  (a)  at  3  hr  with 
contour  Interval  of  O.l,  (b)  at  9  hr  with  contour  interval  of  0.3,  and  (c)  at 
17  hr  with  contour  interval  of  0.3. 
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0 '  subsequently  generates  another  activity  center  of  a  mesolow  which  results  in 
the  w  field  at  11,  13,  and  15  hr,  as  shown  in  Fig.  5.10.  In  Fig.  5.10,  two 
features  are  important:  one  is  the  downward  motion  below  the  upward  motion  area 
in  the  central  portion  of  the  domain  (it  is  especially  evident  in  Fig.  5.10c); 
another  is  the  secondary  maximum  of  upward  motion  beside  the  original  upward 
motion  area.  The  downward  motion  can  reduce  the  heating  ability  of  the  A-S 
scheme  for  the  upper  atmosphere,  as  illustrated  in  Fig.  5.11a,  which  shows  the 
profiles  of  the  adiabatic  cooling  versus  convective  heating  for  the  central 
column  at  13  hr.  It  can  be  seen  that  the  cooling  rate  exceeds  the  heating  rate 
for  the  atmosphere  above  350  mb  so  that  the  temperature  would  decrease 
afterwards  above  that  level.  This  distribution  is  attributed  to  the  total 
vertical  velocity  shown  in  Fig.  5.11b  where  no  distinct  lower  maximum  is  found 
because  of  the  downward  motion  induced  by  the  model.  Without  sufficient 
low-level  mass  convergence,  the  activities  of  the  high  clouds  are  reduced. 
Another  unfavorable  factor  in  reducing  the  original  convection  around  the 
center  is  Chat  the  secondary  upward  motion  consumes  most  of  the  available 
moisture.  On  the  other  hand,  it  is  found  that  the  secondary  upward  motion  area 
cannot  develop  into  a  well-defined  convective  system  (see  Fig.  5.7c).  This 
nondevelopment  can  also  be  attributed  to  insufficient  low-level  mass 
convergence  because  the  magnitudes  of  the  Imposed  Wg  below  the  secondary  upward 
motion  are  negligibly  small  (due  to  the  property  of  the  sinusoidal  function). 

Fig.  5.12  shows  the  moisture  deviation  from  the  initial  state  (q')  at  the 
different  stages.  The  positive  anomalies  at  3  and  9  hr  are  mostly  due  to  the 
detrainment  of  cloud  air.  The  evaporation  of  detrained  liquid  water  also 
contributes  to  moistening  the  atmosphere.  It  is  seen  in  Fig.  5.12b  that  the 
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Fig.  5.10.  Cross  sections  of  vertical  velocity  (w)  in  units  of  cm  s”^  in  the 
control  run:  (a)  at  11  hr  with  contour  Interval  of  2.0,  (b)  at  13  hr  with 
contour  Interval  of  2.0,  and  (c)  at  15  hr  with  contour  interval  of  1.0. 
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Fig.  5.12.  Cross  sections  of  the  deviation  of  mixing  ratio  (q')  from  the 
initial  condition  in  units  of  g  kg~^  in  the  control  run:  (a)  at  3  hr  with 
contour  Interval  of  0.1,  (b)  at  9  hr  with  contour  Interval  of  0.2,  and  (c)  at 
17  hr  with  contour  interval  of  0.2. 
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net  Increase  In  the  mixing  ratio  after  the  9-hour  simulation  is  rather  small 
with  the  maximum  about  1.3  g  Kg“^.  The  negative  anomalies  outside  the 
convective  area  are  undoubtedly  due  to  the  downward  motion •  The  decrease  in 
the  positive  anomalies  around  the  lower  center  from  9  to  17  hr  is  also  due  to 
the  downward  motion.  Fig.  3.13  shows  the  distributions  of  the  cloud  base  mass 
fluxes  (mg)  as  functions  of  cloud  top  height  and  distance  from  the  center  at 
the  different  stages.  At  3  hr,  blmodal  distributions  are  clearly  present  in 
the  convective  area.  However,  in  the  mature  stage,  only  intense  high  clouds 
are  found  in  the  central  portion  of  the  domain.  Finally,  in  the  decaying 
stage,  the  cloud  population  becomes  much  less. 

The  evolution  of  the  mesoscale  fields  in  the  period  of  24  hours  has  been 
discussed  in  this  section.  The  life  cycle  of  this  simulated  convective  system 
may  be  defined  by  categorizing  the  period  from  0  to  5  hr  as  the  developing 
stage,  the  one  from  6  to  12  hr  as  the  mature  stage,  and  the  one  from  13  to  19 
hr  as  the  decaying  stage.  The  period  from  20  to  24  hr  is  practically  in  a 
steady  state.  In  the  next  section,  we  shall  compare  the  model  results  with  the 
observational  studies  of  this  cloud  band  event. 

5.1.3  Comparison  with  the  observations: 

Figure  5.14  adapted  from  Ogura  et  al.  (1979)  shows  the  vertical  velocity 
(u))  at  the  center  of  the  cloud  band  as  a  function  of  time  and  height,  which  is 
determined  by  the  wind  data  from  the  observational  network  with  a  few  hundred 
kilometers  of  station  separation.  The  main  features  of  this  ui  field  are 
described  as  follows.  Prior  to  and  during  the  initial  stage,  low-level  upward 
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Fig.  5.13.  Cross  sections  of  the  cloud  base  mass  flux  (mR)  In  units  In  mb  hr~^ 
In  the  control  run:  (a)  at  3  hr  with  contour  interval  of  0.8,  (b)  at  9  hr  with 
contour  interval  of  1,0,  and  (c)  at  17  hr  with  contour  interval  of  0.4. 
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Fig.  5.13.  (Continued) 
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motlon  was  present  while  weak  subsidence  was  present  at  middle  and  upper 
levels.  As  convection  developed,  u>  Increased  In  Intensity  and  upward  motion 
was  present  throughout  a  deep  layer.  The  local  maximum  was  located  around  the 
400  mb  level  during  the  mature  stage  of  the  cloud  band  at  15  GMT.  The  decaying 
stage  was  characterized  by  rapid  weakening  of  the  ascending  motion  at  low 
levels  followed  by  the  development  of  descending  motion,  while  the  ascending 
motion  with  a  substantial  Intensity  remained  at  upper  levels.  Figure  5.15 
shows  the  evolution  of  the  model-predicted  w  field  horizontally  averaged  over 
160  km  In  the  central  portion  of  the  domain,  which  is  considered  the  width  of 
the  cloud  band  simulated  by  the  model.  It  can  be  seen  that  the  model  result 
generally  agrees  well  with  the  observation  shown  in  Fig.  5.14.  The  area 
averaged  w  (w)  shown  In  Fig.  5.15  also  increases  In  intensity  with  time  and 
reaches  the  maximum  of  about  15  cm  s”^  at  9  hr,  which  is  somewhat  smaller  than 
the  maximum  shown  In  Fig.  5.14.  The  observation  indicates  that  the  height  of 
the  local  maximum  of  m  Increases  with  time  from  0  to  9  GMT.  If  we  define  this 
9  hour  period  to  be  the  developing  stage  of  the  cloud  band,  then  that  of  the 
model  results  only  takes  4  hours.  As  discussed  previously,  this  Is  because  the 
A-S  scheme  generates  high  clouds  much  faster  than  the  real  atmosphere.  For  the 
decaying  stage,  the  characteristics  of  the  low-level  downward  motion  are  also 
clearly  observed  in  Fig.  5.15. 

Figure  5.16  shows  the  model-predicted  surface  precipitation  as  a  function 
of  time  and  distance  from  the  center.  The  width  and  distribution  of 
significant  precipitation  almost  coincide  with  those  of  upward  motion  as  they 
should  be.  The  maximum  intensity,  about  3  mm  hr”^ ,  is  located  at  the  center  of 
the  domain  at  8  hr.  Figure  5.17  shows  the  observed  precipitation  accumulation 
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Fig.  5.16.  Model-predicted  rainfall  accumulation  in  the  control 
function  of  time  and  distance  from  the  center  In  units  of  mm  hr~^  with 
interval  of  0.3. 
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for  the  cloud  band  from  12  to  18  GMT  12  August  1974.  It  is  found  that  the 
model  result  agrees  well  with  the  observation  both  in  width  scale  and  intensity 
of  the  precipitation. 


The  observational  studies  indicate  that  both  the  net  heating  rate  by  a 

well-developed  cloud  cluster  and  the  adiabatic  cooling  rate  by  the  large-scale 

lifting  are  typically  10  K  day“^  in  order  of  magnitude.  In  other  words,  the 

d  6 

temperature  change  ( — .]  in  the  tropical  region  is  a  small  residue  between  the 
two  large  terms  of  cooling  by  the  large-scale  lifting  and  heating  by  the  clouds 
and  is  1  K  day”'^  in  order  of  magnitude.  The  model  presented  here  appears  to  be 
capable  of  simulating  this  very  unique  feature  in  the  tropics  as  Illustrated  in 
Table  5.1  which  presents  the  deviation  of  area-averaged  potential  temperature 
and  mixing  ratio  at  9  hr  from  the  initial  states.  This  implies  that  the  A-S 
scheme  does  have  a  potential  usefulness  in  predicting  tropical  phenomena. 


5.2  Sensitivity  tests: 


5.2.1  The  deletion  of  the  Imposed  large-scale  lifting: 


The  model  results  presented  In  the  preceding  section  confirm  the 
observational  fact  that  the  low-level  large-scale  lifting  is  necessary  to 
trigger  the  formation  of  an  organized  convective  system  even  if  the  atmospheric 
stratification  is  conditionally  unstable.  Here  we  shall  undertake  an 


experiment  by  deleting  the  Imposed  low-level  lifting  after  9  hours  of 
simulation  time  to  investigate  the  possible  situation  that  the  large-scale 
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Height  (Km] 

7’  (k) 

T  (g  Kg  ) 

0.45 

0.04 

0.15 

1.35 

-0.22 

0.42 

2.25 

-0.16 

0.24 

3.15 

0.21 

0.15 

4.05 

0.75 

0.25 

4  .95 

0.66 

0.73 

5  .85 

0  .61 

0.03 

6  .75 

0.64 

-0.06 

7.65 

0  .95 

-0.02 

8.55 

1  .15 

-0.05 

9.45 

1.14 

-0.10 

10.35 

1.38 

-0.07 

11.25 

1.36 

-0.03 

12.15 

0.82 

0.07 

13.05 

0  .07 

0.05 

13.95 

0.09 

0.01 

Table  5.1.  The  deviation  of  area-averaged  potential  temperature  and  mixing 
ratio  from  the  Initial  conditions  In  the  control  run. 


lifting  goes  out  of  the  convective  area  (or  dies  out)  after  the  development  of 
the  cloud  cluster. 

Figure  5.18  shows  the  w  fields  at  11,  13,  and  15  hr  after  the  deletion  of 
Wq  at  9  hr.  The  result  clearly  indicates  a  drastic  decrease  in  the  maximum  of 
w  from  17  cm  s”^  at  9  hr  (see  Fig.  5.7b)  to  3  cm  s“^  at  15  hr.  It  is  also 
clearly  observed  from  Fig.  5.18  that  the  original  convective  area  is  cut  into 
halves  by  the  intrusion  of  downward  motion  from  the  low  levels.  One  of  them 
rapidly  decays  and  the  other  propagates  laterally.  Figure  5.19  shows  6'  fields 
at  11,  13,  and  15  hr.  The  magnitudes  of  the  positive  anomalies  around  the 
center  are  smaller  while  new  maxima  are  generated  far  away  from  the  center. 
The  original  negative  anomaly  at  the  low  levels  around  the  center,  as  shown  in 
Fig.  5.8b,  disappears  because  of  the  deletion  of  Wq  and  its  associated 
adiabatic  cooling.  The  negative  anomaly  at  the  upper  levels  around  the  center, 
shown  in  Fig.  5.19,  is  due  to  the  smaller  heating  ability  of  the  A-S  scheme, 
since  the  deletion  of  Wq  decreases  the  population  of  the  high  clouds.  In  fact, 
the  basic  reason  for  initiating  the  rapid  dissipitation  in  this  run  can  be 
explained  by  Fig.  5.20  which  bears  a  similar  rationale  with  that  of  Fig.  5.11 
for  the  decaying  stage  in  the  control  run.  The  deletion  of  the  Wq  makes  the 
lower  maximum  of  the  total  vertical  velocity  disappear,  as  shown  in  Fig.  5.20b, 
and  the  adiabatic  cooling  rate  exceed  the  convective  heating  rate  for  the 
middle  and  upper  atmosphere,  as  shown  in  Fig.  5.20a.  Meanwhile,  since  the 
adiabatic  cooling  rate  due  to  Wq  at  the  low  levels  does  not  exist  any  longer,  a 
large  excess  of  convective  heating  over  adiabatic  cooling  at  the  low  levels  is 
observed.  Obviously,  this  Implies  that  the  atmospheric  static  stability  is 
destabilized.  Nonetheless,  from  the  viewpoint  of  the  coupling  or  Interactions 
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Fig.  5.18.  Cross  sections  of  vertical  velocity  (w)  in  units  of  cm  in  Run 

2:  (a)  at  11  hr  with  contour  Interval  of  1.0,  (b)  at  13  hr  with  contour 
interval  of  1.0,  and  (c)  at  15  hr  with  contour  interval  of  0.7. 
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Fig.  5.19.  Cross  sections  of  the  deviation  of  potential  temperature  (0' 
the  Initial  condition  In  units  of  K  in  Run  2:  (a)  at  11  hr  with 
Interval  of  0.2,  (b)  at  13  hr  with  contour  interval  of  0.2,  and  (c)  at 
with  contour  Interval  of  0.2. 
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Fig.  5.19.  (Continued) 


between  Che  mesoscale  circulation  and  cumulus  convection,  the  original 
horizontal  temperature  gradient  around  Che  center  is  reduced  because  of  the 
distributions  of  heating  and  cooling,  shown  in  Fig.  5.20a,  and  its 
corresponding  upward  motion  is  consequently  weakened.  It  should  be  noted  that 
several  aspects  of  the  decaying  stages  of  the  two  runs  are  different.  Other 
than  the  difference  of  the  time  scales  of  dissipation,  there  is  no  distinct 
propagation  of  w  fields  in  the  control  run,  as  shown  in  Fig.  5.10,  while  rapid 
propagation  is  evident  in  Fig.  5.18.  This  is  because  in  Che  control  run  the 
imposed  Wq  Is  always  present .  Although  the  imposed  Wq  Is  partially  compensated 
by  the  downward  motion  in  the  decaying  stage,  the  net  effect  of  Wq  is  still  to 
release  Che  potential  Instability  of  Che  lower  atmosphere  and  practically 
serves  as  a  continuous  heating  source  through  the  A-S  scheme.  It  is  believed 
that  suddenly  deleting  this  heating  source  can  result  in  the  rapid  lateral 
propagation  of  the  convective  system.  More  relevant  discussion  along  the 
problem  of  system  propagation  will  be  given  in  Section  5.2.4  where  Kuo's  scheme 
is  Incorporated  into  the  model.  Figure  5.21  shows  the  vector  fields  at  13  hr 
for  the  two  different  runs.  In  the  control  run,  a  low-level  reverse 
circulation  is  present  with  fairly  strong  inflows  lying  above,  as  shown  in 
Fig.  5.21a,  which  introduces  a  possible  moisture  supply  for  the  secondary 
upward  motion  area.  Obviously,  this  feature  cannot  be  observed  from 
Fig.  5.21b.  Generally  speaking,  the  experiment  presented  in  this  section 
reveals  another  possible  decaying  situation  whose  time  scale  is  much  shorter 
and  in  which  physical  processes  are  much  more  straightforward. 
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Fig.  5.21.  Cross  sections  of  vector  field  of  [u,w]:  (a)  at  13  hr  in  the 
control  run  and  (b)  at  13  hr  In  Run  2.  The  magnitudes  of  w  have  been  enlarged 
by  10  times. 
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5.2.2  The  effects  of  vertical  wind  shear: 


The  initial  horizontal  wind  (Ug)  profile  used  in  this  run  is  shown  in 

Fig.  5.22.  According  to  Eq.  (4.8),  the  corresponding  meridional  temperature 

gradient  f — is  about  -0.4  x  K  km~^  and  this  value  will  be  used  as  a 

3y'' 

constant  throughout  the  integration.  First  of  all,  it  should  be  expected  that 
the  model  results  would  not  stay  symmetric  about  the  center  line  of  the  domain. 
Figure  5.23  shows  the  w  fields  for  the  whole  domain  at  the  different  stages. 
It  is  seen  that  the  convection  does  shift  to  form  in  the  downwind  area  but  not 
very  distant  from  the  center,  with  the  width  of  the  convective  area  about  the 
same  as  that  of  the  control  run  shown  in  Fig.  5.7.  The  main  difference  between 
Figs.  5.23  and  5.7  lies  in  the  distributions  of  downward  motion.  Combining 
Figs.  5.23b  and  5.23c  with  the  intermediate  w  fields  at  11,  13,  and  15  hr  shown 
in  Fig.  5.24,  it  appears  that  the  gradual  intrusion  of  the  downward  motion  from 
the  low  levels  at  9  hr  into  the  right  half  of  the  domain  characterizes  the 
dissipation  of  the  system.  Figures  5.25  and  5.26  show  the  fields  of  9'  and  its 
horizontal  gradient  at  the  different  stages.  Small  negative  temperature 
anomalies  are  found  at  the  top  of  the  convective  portion,  shown  in  Fig.  5.25, 
which  implies  that  the  heating  due  to  high  clouds  is  less  than  that  in  the 
control  simulation.  This  difference  is  caused  by  the  fact  that  the  low-level 


large-scale  lifting  does  not  shift  with  the  horizontal  wind  so  that  the 


magnitudes  of  Wq  under  the  convective  portion  in  the  present  experiment  are 
smaller  than  those  in  the  control  run.  The  temperature  gradient,  shown  in 
Fig.  5.26,  in  conjunction  with  the  ageostrophic  wind  (v^g)  shown,  in  Fig.  5.27, 
explains  the  evolution  of  w  fields  shown  in  Fig.  5.23.  Figure  5.28  shows  the 
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Fig.  5.23.  Cross  Sections  of  vertical  velocity  (w)  in  units  of  cm  in  Run 
3:  (a)  at  3  hr  with  contour  interval  of  2.0,  (b)  at  9  hr  with  contour  interval 
of  2.0,  and  (c)  at  17  hr  with  ocntour  Interval  of  1.0. 
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Fig.  5,23.  (Continued) 
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Fig.  5.24.  Cross  sections  of  vertical  velocity  (w)  in  units  of  cm  s 
3:  (a)  at  11  hr  with  contour  interval  of  2.0,  (b)  at  13  hr  with 
interval  of  2.0,  and  (c)  at  15  hr  with  contour  interval  of  1.0. 
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Fig.  5.25.  Cross  sections  of  the  deviation  of  potential  temperature  (9')  from 
the  initial  condition  in  units  of  K  in  Run  3:  (a)  at  3  hr  with  contour  interval 
of  1.0,  (b)  at  9  hr  with  contour  Interval  of  0.3,  and  (c)  at  17  hr  with  contour 
interval  of  0.3. 
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Fig.  5.25.  (Continued) 
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Fig.  5.26.  Cross  sections  of  the  horizontal  gradient  of  potential  temperature 
in  units  of  K  100  km  in  Run  3:  (a)  at  3  hr  with  contour  interval  of  0.1, 
(b)  at  9  hr  with  contour  Interval  of  0.2,  and  (c)  at  17  hr  with  contour 
Interval  of  0.3. 
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Fig.  5.27.  Cross  sections  of  the  v  ageostrophlc  wind  (v^g)  in  units  of  cm  s~^ 
in  Run  3:  (a)  at  3  hr  with  contour  interval  of  400,  (b)  at  9  hr  with  contour 
interval  of  400,  and  (c)  at  17  hr  with  contour  interval  of  300. 
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Fig.  5.28.  Time-height  section  of  the  model-predicted  vertical  velocity  in  Run 
3  averaged  over  160  km  around  the  center  of  the  model  domain  in  units  of  cm  s~^ 
with  contour  interval  of  2.0. 
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horlzontally  averaged  w  (w)  over  160  km  around  the  central  portion  of  the 
domain,  which  covers  40  km  In  the  left  half  and  120  km  in  the  right  half.  No 
substantial  differences  are  found  In  the  life  cycle  and  Intensity  of  w  when 
compared  with  the  result  of  the  control  run  shown  In  Fig.  5.15.  Generally 
speaking,  the  effects  of  vertical  wind  shear  may  be  of  no  Importance  In  this 
model  framework. 

Many  cloud  modeling  studies  (Schlesslnger ,  1973;  Wllhelmson,  1974;  Thorpe 
et  al. ,  1982;  and  many  others)  have  concluded  that  the  vertical  wind  shear  is 
responsible  for  the  long-lived,  multi-cell  thunderstorms.  The  shortness  of  the 
life  span  of  individual  thunderstorms  results  from  the  process  that  the 
downdrafts,  which  develop  Inside  the  cell  due  to  water  loading,  cut  off  the 
supporting  upward  motion.  In  contrast,  the  vertically-sheared  flow  can 
effectively  dislocate  the  water  particles  from  the  updraft  area.  This  effect 
makes  the  thunderstorms  live  longer  or  turn  Into  the  multi-cell  storms.  With 
this  recognition,  the  insensitivity  of  this  model  to  the  vertical  wind  shear  is 
not  surprising,  since  no  precipitating  particles  are  Included  in  the  present 
model . 


5.2.3  The  effects  of  cumulus  downdrafts: 


The  semi -prognostic  studies  with  the  Inclusion  of  cumulus  downdrafts  have 
been  presented  in  Section  3.4.  The  results  Indicated  that  the  agreement 
between  the  predicted  and  observed  cloud  heating  effects  is  as  good  as  that  In 
the  case  without  the  downdrafts.  Nevertheless,  it  is  not  guaranteed  that  the 
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results  of  the  mesoscale  model  would  be  insensitive  to  the  downdraft  effects 
because  of  the  following  reasons.  First,  from  the  discussion  in  Section  5.1, 
it  is  recognized  that  In  this  mesoscale  model  the  development  of  the  convective 
system  is  exclusively  initiated  by  the  small  difference  between  the  convective 
heating  and  adiabatic  cooling.  The  inclusion  of  the  cumulus  downdraft  effects 
may  significantly  alter  this  "small  difference"  and  consequently  result  in  very 
different  mesoscale  fields.  Secondly,  the  feedback  mechanism  between  the  cloud 
layer  and  subcloud  layer  through  the  cumulus  downdrafts  may  play  an  important 
role  in  the  simulation.  As  discussed  earlier,  the  detralnment  of  downdrafts 
into  the  mixed  layer  and  the  evaporation  of  the  precipitation  in  the  mixed 
layer  can  Introduce  the  modification  of  the  mixed  layer  thermodynamic 

properties  (9^^  and  q^) ,  which  are  important  factors  in  determining  the  cloud 
properties.  Here  it  is  worthwhile  to  mention  that,  in  the  control  run,  the 
mixed  layer  properties  are  only  governed  by  the  large-scale  horizontal 
advection  and  turbulent  eddy  fluxes  so  that  the  time  rate  of  change  of  9j^  is  as 
small  as  10“^  K  day“^  in  order  of  magnitude. 

It  should  be  mentioned  that  e  =  -0.3  (see  Eq.  (3.19))  is  used  in  the 

present  experiment.  Figure  5.29  shows  the  potential  temperature  deviation  (0') 

from  the  initial  state  at  3,  6,  9,  and  12  hr.  As  compared  with  Fig.  5.8,  it  is 
seen  that  the  maximum  temperature  drop  in  the  mixed  layer  is  as  large  as  0.7  K 
at  3  hr  due  to  the  downdrafts.  It  is  also  found  that  a  considerable 
I  temperature  decrease  is  present  in  the  low  levels  around  the  center  above  the 

mixed  layer,  which  can  be  expected  to  Induce  the  low-level  reverse  circulation 
more  quickly  and  strongly  than  in  the  control  run.  Furthermore,  the  maximum 
I  temperature  increase  due  to  convective  heating  is  only  about  0.9  K  at  9  hr. 
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Fig.  5.29.  Cross  sections  of  the  deviation  of  potential  temperature  (0')  from 
the  initial  condition  in  units  of  K  in  Run  U:  (a)  at  3  hr  with  contour  interval 
of  O.l,  (b)  at  6  hr  with  contour  Interval  of  0.2,  (c)  at  9  hr  with  contour 
Interval  of  0.2,  and  (d)  at  12  hr  with  contour  Interval  of  0.3. 
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The  basic  reason  for  Che  situation  described  above  can  be  demonstrated  in 
Fig.  5.30,  which  is  the  counterpart  of  Fig.  5.5,  showing  the  difference  in 
magnitudes  between  the  convective  heating  and  adiabatic  cooling  for  the  very 
beginning  of  the  simulation.  It  is  found  In  Fig.  5.30  that  the  excess  of  the 
adiabatic  cooling  over  the  convective  heating  in  the  lower  atmosphere  is 
slightly  larger  than  that  in  Fig.  5.5  while  the  excess  of  the  convective 
heating  over  the  adiabatic  cooling  in  the  upper  atmosphere  is  slightly  smaller. 
Obviously,  these  differences  are  attributed  to  the  modification  of  the  A-S 
scheme  due  to  the  Inclusion  of  the  downdraft  effects  In  the  cloud  layer.  Based 
on  the  discussion  of  Figs.  5.5  and  5.30,  two  important  lessons  are  learned. 
First,  in  the  present  mesoscale  model,  the  early  development  of  the  convective 
system  greatly  depends  upon  the  difference  between  the  magnitudes  of  the 
convective  heating  and  adiabatic  cooling.  Secondly,  we  may  further  confirm  our 
previous  statement  about  the  different  purposes  between  the  semi-prognostic 
approach  and  fully  prognostic  model  in  testing  a  cumulus  parameterization 
scheme.  The  small  difference  between  the  two  profiles  in  Fig.  5.30  or  Fig.  5.5 
may  be  of  no  concern  in  the  semi-prognostic  approach  but  can  be  extremely 
crucial  in  a  fully  prognostic  model,  especially  in  simulating  the  development 
of  a  mesoscale  convective  system. 

It  is  obvious  that  the  weaker  convective  activities  in  this  run  are  also 
attributed  to  the  temperature  drop  of  the  mixed  layer  resulting  from  the 
downdraft  effects  in  terms  of  the  downdraft  heat  flux  and  the  evaporation  of 
precipitation  (see  Eq.  (4.28)).  The  result  shows  that  the  contributions  of 
these  two  factors  to  the  decrease  of  are  about  the  same  in  magnitude.  As 
indicated  by  Frank  (1980)  in  an  observational  study  of  GATE  convective  systems. 
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Fig.  5.30.  Profiles  of  convective  heating  (dashed  line)  in  response  to  the 
imposed  Wq  and  adiabatic  cooling  (solid  line)  due  to  the  Imposed  Wq  in  Run  4. 
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when  the  boundary  layer  is  disturbed  by  downdrafts,  the  convective  heating  in 
the  middle  and  upper  troposphere  is  generally  weaker  than  the  boundary  layer 
cooling  and  plays  a  minor  role  in  stabilizing  the  atmosphere.  Observations 
also  indicate  that  the  stabilization  resulting  from  downdrafts  is  sufficient  to 
inhibit  further  convective  development  until  the  boundary  layer  recovers 
(Fitzjarrald  and  Garstang,  1981).  In  the  preceding  paragraph,  we  have  shown 
the  effect  of  the  modification  of  the  A-S  scheme  due  to  the  inclusion  of 
cumulus  downdrafts  on  the  heating  profile.  Here  we  shall  illustrate  the 
sensitivity  of  the  heating  profile  to  the  temperature  in  the  mixed  layer.  To 
do  so,  we  select  the  center  column  at  3  hr  and  compare  its  heating  profile  with 
the  one  which  is  computed  from  the  same  atmospheric  conditions  except  for  the 
mixed  layer  temperature,  which  is  1  K  higher.  The  result  shown  in  Fig.  5.31 
Indicates  that  a  significant  difference  is  apparently  present  between  750  to 
550  mb  and  below  850  mb  and  that  the  case  with  thr  lower  mixed  layer 
temperature  provides  less  convective  heating.  With  the  recognition  that  the 
model  results  are  sensitive  to  the  difference  between  the  convective  heating 
and  adiabatic  cooling  at  the  initial  stage,  as  discussed  earlier,  another 
implication  of  the  result  shown  in  Fig.  5.31  is  that  the  evolution  of  the 
mesoscale  fields  is  sensitive  to  the  initial  mixed  layer  thermodynamic 
properties . 

Figure  5.32  shows  the  corresponding  w  fields  at  3,  6,  9,  and  12  hr.  As 
compared  with  Fig.  5.7,  significant  differences  are  found  in  that  (a)  the 
magnitudes  of  w  are  smaller  and  the  maximum  of  w  in  the  entire  simulation  is 
only  about  10  cm  s”^  at  9  hr,  (b)  the  upward  motion  area  slightly  shifts 
outward  with  time,  and  (c)  Intense  downward  motion  occupies  all  of  the  central 
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Fig.  5.32.  Cross  sections  of  the  vertical  velocity  (w)  In  units  of  cm  s~^  in 
Run  4:  (a)  at  3  hr  with  contour  Interval  of  0.6,  (b)  at  6  hr  with  contour 
Interval  of  l.O,  (c)  at  9  hr  with  contour  Interval  of  1.0,  and  (d)  at  12  hr 
with  contour  interval  of  1.0. 


columns  at  9  and  12  hr.  In  passing.  It  should  be  mentioned  that  no  distinct 
change  In  the  mixing  ratio  In  the  mixed  layer  (<!]«()  Is  observed  because  as  the 
downdrafts  bring  the  drier  air  Into  the  mixed  layer,  the  precipitation 
evaporates  to  moisten  the  mixed  layer.  Evidently,  the  experiment  with  the 
cumulus  downdraft  effects  presented  above  has  provided  less  realistic  results 
than  the  case  without  the  downdrafts.  Nevertheless,  serving  as  a  conceptional 
experiment,  this  simulation  does  reveal  several  aspects  of  the  cumulus 
downdraft  effects  not  only  In  conjunction  with  the  AS  scheme  but  also  with  the 
simulation  of  a  mesoacale  convective  system. 

5.2.4  The  teat  of  Kuo's  scheme: 

Soong  and  Tao  (1980)  made  a  comparison  between  the  results  from  explicit 
calculations  for  clouds  and  that  derived  from  Kuo's  1965  cumulus 
parameterization  scheme  within  the  framework  of  the  cumulus  ensemble  model 
discussed  previously.  They  found  that  the  cloud  heating  profile  computed  by 
Kuo's  scheme  had  a  maximum  In  the  upper  atmosphere  while  the  heating  rate  below 
700  mb  was  negligibly  small,  which  does  not  agree  with  the  observed  heating 
profile  at  all.  In  the  present  experiment,  we  shall  replace  the  A-S  scheme  by 
Kuo's  1974  scheme  (Kuo,  1974)  with  the  parameter  b  set  to  0,  as  suggested  by 
Kuo  (1974)  for  highly  convective  regions.  Setting  b  to  0  Implies  that  all  of 
the  moisture  convergence  into  an  atmospheric  column  will  precipitate  to  heat 
the  atmosphere.  It  would  be  interesting  to  see  whether  or  not,  or  to  what 
extent,  Kuo's  scheme  is  able  to  simulate  the  development  of  a  mesoscale 
convective  system  In  the  present  model  framework. 
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The  detailed  derivation  of  Kuo's  scheme  has  been  given  In  Kuo  (1965, 
1974).  Here,  we  briefly  summarize  the  scheme  as  follows.  The  total  rate  of 
moisture  accession  per  unit  horizontal  area  (M^)  is  given  by: 


Zrr 

M,.  =  /  V»(qV)  pdz  +  surface  evaporation  rate 
0 


(5.3) 


where  Zj  is  the  top  of  the  domain  and  the  surface  evaporation  rate  can  be 
expressed  by  (4.37)  or  another  bulk  formula.  Having  assumed  that  all  of  the 
moisture  convergence  Is  for  heating  the  atmosphere,  we  can  determine  the  rate 
of  cloud  production  per  unit  time  (Q^)  as: 


%  - 


where  is  expressed  by: 


(5.4) 


(5.5) 


where  z^^  represents  the  height  of  cloud  top.  Equation  (5.5)  represents  the 
amount  of  moisture  which  will  condense  and  precipitate  to  raise  the 
environmental  temperature  (T^)  to  cloud  temperature  (T^).  T^.  has  the  moist 
adiabatic  value  and  can  be  obtained  by  the  saturation  technique.  Now  assume 
that  the  prediction  is  made  for  a  time-step  At  giving  a  temperature  T*  without 
the  effects  of  cumuli  at  t  +  At  as: 


T(t  -  At)  +  2At*Fj. 


-  •  •  •  •  •*-  •  •*'.  •*„  w*. 


(5.6) 


where  represents  the  advectlon  terms  in  the  thermodynamic  equation.  Using 
T  and  Its  corresponding  to  obtain  and  using  the  current  V«qV  and  surface 
evaporation  to  obtain  we  can  determine  by  (5.4).  Then,  an  adjusted 
temperature  T(t  +  At)  Is  determined  by: 

T(t  +  At)  -  T*  +  At(T^  -  T*)  (5.7) 

Note  that,  since  the  value  of  Is  at  time  t,  only  one  At  Is  multiplied  in 
(5.7).  Accordingly,  the  rate  of  precipitation  per  unit  time  per  unit  volume  in 
the  column  Is  Q^.  Cp  (T^.  -  T*)/L.  Since  b  Is  set  to  0,  there  is  no  Increase  or 
decrease  of  humidity  as  the  convection  takes  place.  During  the  Integration  of 
the  model.  If  the  computed  amount  of  exceeds  the  specified  value  of 
0.3  X  10**^,  the  convection  Is  assumed  to  take  place. 

It  should  be  mentioned  that  in  this  model  in  (5.7)  Is  determined  from 
the  non-entraining  plume  such  that  only  the  highest  clouds  are  allowed  to  form. 
Figure  5.33  shows  the  evolution  of  w  fields  for  the  period  from  1  to  5  hr  at 
every  1  hr  Interval.  At  3  hr,  the  Intensity  of  w  reaches  its  maximum  with  the 
magnitude  about  10  cm  s~^,  which  Is  considered  significantly  smaller  In 
comparison  with  the  observation  and  the  result  simulated  In  the  control  run. 
The  corresponding  9'  fields  are  shown  In  Fig.  5.34.  The  most  distinct  feature 
In  Figs.  5.33  and  5.34  Is  the  rapid  lateral  propagation  of  the  convective 
system  with  the  speed  about  30  m  s“^.  As  compared  with  the  results  In  the 
control  run.  It  Is  felt  that  the  fast  lateral  propagation  In  the  present 
experiment  Is  mainly  due  to  the  heating  profile  generated  by  Kuo's  scheme.  As 
Indicated  previously,  in  the  developing  stage  of  the  control  run,  the  vertical 
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Fig.  5.34.  Cross  sections  of  the  deviation  of  potential  temperature  (0')  in 
units  of  K  in  Run  5  using  Kuo's  scheme:  (a)  at  1  hr  with  contour  interval  of 
O.l,  (b)  at  2  hr  with  contour  interval  of  0.1,  (c)  at  3  hr  with  contour 
interval  of  0.2,  (d)  at  4  hr  with  contour  Interval  of  0.2,  and  (e)  at  5  hr  with 
contour  Interval  of  0.2. 
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Fig.  5.34.  (Continued) 
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heating  profiles  have  Che  double  maximum  distribution,  as  shovn  in  Fig.  5.9a. 
Figure  5.35  shows  the  vertical  profile  of  the  convective  heating  rate  for  the 
central  column  at  1  hr.  From  Eq.  (5.7),  It  is  seen  that  Che  amount  of 
convective  heating  Is  proportional  to  (T^  -  T*)  which  ordinarily  has  the 
maximum  at  upper  levels,  as  illustrated  by  Fig.  5.35.  On  the  other  hand,  the 
negligibly  small  heating  in  the  lower  atmosphere  cannot  compensate  the 
adiabatic  cooling  due  Co  the  imposed  vertical  motion  (wq) .  Accordingly,  this 
results  in  the  large  amount  of  negative  temperature  anomalies  around  Che  lower 
central  portion,  as  shown  in  Fig.  5.34. 

It  is  generally  believed  Chat  a  large  amount  of  upper  level  heating  tends 
Co  excite  gravity  waves  and  causes  the  lateral  propagation.  The  fast  lateral 
dispersal  of  the  convective  heating  has  been  discussed  in  a  review  paper  by 
Frank  (1983).  He  indicated  that  when  the  atmosphere  is  heated  on  a  dynamically 
small  scale  (i.e.,  Che  spatial  horizontal  scale  is  smaller  than  the  Rossby 
radius  of  deformation),  such  as  Che  release  of  latent  heat  in  a  convective 
mesoscale  system,  the  divergence  field  is  usually  not  significantly  constrained 
due  to  Che  lack  of  stiffness  of  rotational  flow  resulting  from  the  weakness  of 
inertial  Instability.  Thus,  one  net  result  is  the  lateral  propagation  of  the 
Induced  warming.  Now  recall  that  the  experiment  presented  in  Section  5.2.1 
also  demonstrates  similar  propagation  after  the  low-level  forcing  (Wg)  is 
deleted  (and,  consequently,  so  is  the  low-level  heating),  as  shown  in 
Fig.  5.18,  while  no  such  propagation  is  observed  in  the  control  simulation. 
Thus,  it  is  felt  that  sufficient  low-level  heating  may  be  indispensable  to 
prevent  the  convective  system  from  propagating. 


CHAPTER  VI 


SUMMARY  AND  CONCLUSION 

In  this  report,  we  first  reviewed  the  Arakawa-Schubert  cumulus 
parameterization  scheme  In  a  framework  that  conceptually  divides  the  mutual 
Interaction  of  the  cumulus  convection  and  larger-scale  disturbances  Into  the 
categories  of  the  large-scale  budget  requirements  and  the  quasl-equlllbrlum 
assumption  of  the  cloud  work  function  as  presented  In  Chapter  2.  The  quasl- 
equlllbrlum  assumption  asserts  a  near  balance  between  the  generation  of  moist 
convective  Instability  due  to  large-scale  processes  and  Its  destruction  due  to 
the  clouds.  Since  the  A-S  scheme  bears  considerable  complexity  and,  thereby, 
different  Implementations  may  produce  quite  different  results,  we  have  been 
cautious  In  selecting  the  computational  algorithm,  as  reported  In  Section  3.2. 
The  algorithm  employed  In  this  report  Is  different  from  that  used  by  Lord 
(1978,  1982)  and  appears  to  be  capable  of  providing  more  accurate  results.  We 
also  reviewed  the  technical  aspects  involved  In  solving  the  cloud  base  mass 
flux  distribution  equation,  which  has  been  restated  by  a  linear  programming 
problem.  A  software  package  was  found  to  be  very  efficient  in  solving  the 
problem. 

We  first  applied  the  A-S  scheme  through  the  semi-prognostic  approach  to 
two  different  data  sets:  one  for  the  tropical  cloud  band  event  and  the  other 
for  the  composite  data  for  GATE  Phases  II  and  III.  Basically,  the  semi¬ 
prognostic  approach  Is  designed  to  Investigate  the  capability  of  a  cumulus 
parameterization  scheme  In  reproducing  the  observed  cloud  heating  and  drying 


profiles.  The  performance  of  the  A-S  scheme  presented  in  Section  3.3  appears 
quite  successful  In  serving  the  purpose.  More  Importantly,  the  results  confirm 
the  potential  advantages  of  the  quasl-equlllbrlum  assumption,  which  Is 
considered  the  core  of  the  A-S  scheme.  Nevertheless,  a  deficiency  of  the  A-S 
scheme  was  found  to  lie  In  the  underestimation  of  the  condensation  and 
evaporation  rates  when  the  semi-prognostic  results  were  compared  with  those 
derived  from  the  cloud  ensemble  model.  Even  though,  from  a  practical  point  of 
view,  this  deficiency  can  be  considered  minor,  the  results  and  the  discussion 
of  the  cumulus  downdraft  effects  In  Section  3.4  have  pointed  out  a  way  to 
Improve  the  A-S  scheme.  In  doing  so,  we  have  modified  the  A-S  scheme  by 
incorporating  the  downdraft  model  formulated  by  Johnson  (1976)  but  with  a 
different  parameter  e  which  is  a  measure  of  the  downdraft  intensity  relative  to 
the  updraft  intensity.  The  results  showed  that  the  modified  A-S  scheme  not 
only  substantially  improved  the  condensation  and  evap  'ration  rates  but  also 
predicted  a  downdraft  mass  flux  which  agreed  fairly  well  with  the  result  of  the 
cloud  ensemble  model.  A  drawback  introduced  by  this  downdraft  model  lies  In 
the  overestimation  of  the  moistening  effect  resulting  from  the  detralnment  of 
the  cloud  air  at  the  cloud  tops.  Because  of  the  lack  of  a  sound  scientific 
basis  to  determine  e ,  it  is  felt  that  this  downdraft  study  through  the  semi¬ 
prognostic  approach  Is  only  an  initial  step  toward  a  more  complete  examination 
of  the  role  played  by  the  cumulus  downdrafts  in  the  cumulus  parameterization 
problem  for  prognostic  purposes. 

A  more  severe  test  of  the  A-S  scheme  was  subsequently  made  by 
incorporating  the  scheme  into  a  two-dimensional  primitive  mesoscale  model  to 
simulate  the  tropical  cloud  band  event  that  was  considered  in  the  previous 


paragraph.  The  details  of  Che  mesoscale  model  have  been  presented  in  Chapter 
4.  The  principal  goals  have  been  to  Investigate  (a)  the  development  of  the 
mesoscale  convective  system  resulting  from  the  convective  heating  as  estimated 
by  the  A-S  scheme,  (b)  the  evolution  of  the  mesoscale  fields  and  its  associated 
physical  processes,  (c)  the  sensitivity  of  the  model  to  different  environmental 
conditions,  and  (d)  the  response  of  the  A-S  scheme  to  the  mesoscale 
circulation. 

Section  5.1  presented  the  results  of  a  24-hour  Integration  for  the  control 
simulation  in  which  the  effects  of  cumulus  downdrafts  and  the  vertical  shear  of 
Che  ambient  wind  were  not  Included.  The  evolution  of  Che  mesoscale  circulation 
in  conjunction  with  the  convective  heating  due  to  cumulus  convection,  the 
adiabatic  cooling  due  to  total  upward  vertical  motion,  and  the  adiabatic 
warming  due  to  mesoscale  downward  motion  was  examined.  Initial  conditions  were 
calm  except  for  the  imposed  large-scale,  low-level  circulation  [wq,  uq]  that 
releases  the  potential  Instability  through  the  applications  of  the  A-S  scheme. 
The  first  6-hour  period  was  dominated  by  the  convective  heating  located  in  the 
central  portion  of  the  domain  where  the  imposed  Wq  was  located.  The  second 
6-hour  period  was  still  basically  dominated  by  the  convective  heating.  During 
this  period,  however,  the  accumulation  of  the  adiabatic  warming  outside  the 
cloud  cluster  area  gradually  became  appreciable.  This  subsequently  generated  a 
reverse  circulation  at  the  low  levels  and  the  decaying  stage  set  in  for  the 
third  6-hour  period.  The  last  6-hour  period  was  more  or  less  in  a  steady 
state.  One  gratifying  aspect  is  that  the  model  results  bear  considerable 
resemblance  to  many  aspects  of  the  observations  including  (a)  the  evolution  of 
area-averaged  vertical  velocity  (w) ,  (b)  the  intensity  and  distribution  of 


surface  precipitation,  (c)  the  width  of  the  mesoscale  convective  system,  and 
(d)  the  development  of  low-level  downward  motion  during  the  decaying  stage. 


Although  the  model  can  provide  the  characteristics  of  the  decay  of  the 
system  in  terms  of  the  rapid  weakening  of  upward  motion  and  the  development  of 
downward  motion  at  the  lower  atmosphere  with  self-consistent  physical 
processes,  the  realism  of  this  decaying  situation  may  be  uncertain  when  one 
recalls  the  following  fact.  It  Is  commonly  documented  that  the  decay  of 
individual  thunderstorms  Is  caused  by  the  low-level  water  particles  which 
reduce  the  buoyancy  force  around  the  cloud  base  and  consequently  Induce 
precipitating  downdrafts  to  cool  and  dry  the  sub-cloud  layer.  Obviously,  all 
of  these  effects  are  not  Included  In  this  particular  simulation.  Nevertheless, 
the  statement  described  above  is  mostly  deduced  from  the  cloud  simulation  with 
explicit  calculation  for  the  cloud-scale  properties  for  an  individual 
convective  entity.  To  date,  the  complete  physical  reasons  for  the  dissipation 
of  the  mesoscale  convective  systems  are  not  well  understood.  If  we  take  a  more 
macroscopic  point  of  view,  the  physical  processes  of  the  decaying  stage  deduced 
from  our  mesoscale  model,  as  discussed  in  Sections  5.1  and  5.2.1,  may  be  valid 
to  some  extent.  This  argument  appears  especially  significant  when  the  decaying 
processes  presented  in  Section  5.2.1  for  run  2  are  considered,  where  we  dealt 
with  a  decaying  situation  by  suddenly  deleting  the  large-scale,  low-level 
forcing. 


A  test  on  the  effects  of  the  vertical  wind  shear  was  presented  in  Section 
5.2.2.  Since  the  water  loading  Is  not  included  in  this  mesoscale  model,  no 
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substantlal  differences  from  Che  control  run  are  observed  except  for  the  slight 
shift  of  the  convective  activities  to  the  downwind  area. 

A  test  on  the  effects  of  Che  cumulus  downdrafts  was  presented  in  Section 
5.2.3.  It  was  found  that  the  A-S  scheme  with  the  incorporation  of  the 
downdrafts  tended  to  stabilize  Che  atmosphere  by  the  relatively  large  excess  of 
the  adiabatic  cooling  over  the  convective  heating  in  the  lower  atmosphere  at 
Che  initial  stage.  In  addition,  the  cooling  of  the  mixed  layer  due  Co  the 
detralnment  of  downdrafts  and  Che  evaporation  of  precipitation  also  reduced  Che 
convective  heating  significantly.  Thus,  less  realistic  results  were  obtained 
when  compared  with  Che  control  simulation.  This  experiment  also  revealed  that 
the  evolution  of  the  mesoscale  fields  is  highly  contingent  upon  the  initial 
mixed  layer  thermodynamic  properties  and  the  difference  between  the  magnitudes 
of  Che  convective  heating  and  adiabatic  cooling  at  Che  beginning  of  the 
simulation.  However,  this  sensitivity  by  no  means  implies  that  Che  mesoscale 
model  presented  in  this  dissertation  is  a  susceptible  one.  As  indicated  by 
Frank  (1983),  in  "dynamically  small"  circulations  such  as  a  mesoscale 
convective  system,  the  timing  and  magnitude  of  Che  latent  heat  release  are 
dominant  factors  in  determining  the  life  cycle  and  structure  of  the  convective 
system  because  Che  divergent  wind  field  is  relatively  strong  in  such  a  system 
and  Is  very  sensitive  to  the  convective  heating. 

An  experiment  was  carried  out  in  Section  5.2.4  by  replacing  the  A-S  scheme 
by  Kuo's  1974  scheme.  The  results  Indicated  that  Kuo's  scheme  gave  a  much  less 
favorable  comparison,  lying  In  the  rapid  lateral  propagation  of  the  convective 
system.  We  attributed  this  phenomenon  to  the  large  amount  of  convective 


heating  in  the  upper  atmosphere  and  the  negligibly  small  amount  in  the  lower 
atmosphere  since  Kuo's  scheme  assumes  chat  the  vertical  heating  is  proportional 
to  the  temperature  difference  between  the  cloud  and  its  environment.  In 
combination  with  the  effect  of  Che  relatively  strong  adiabatic  cooling  at  Che 
low  levels  due  to  Che  Imposed  Wq,  the  net  result  was  the  rapid  lateral 
dispersal  of  the  convective  system. 

In  essence,  this  study  examined  the  Arakawa-Schubert  cumulus 
parameterization  scheme  through  the  semi-prognostic  approach  and  the  mesoscale 
simulation  of  a  tropical  convective  system.  Both  applications  demonstrated 
quite  encouraging  results.  The  former  served  as  a  basic  test  in  evaluating  the 
scheme.  Without  its  satisfactory  performance,  it  would  have  been  meaningless 
to  further  incorporate  the  scheme  into  any  fully  prognostic  model.  The 
mesoscale  simulation  presented  here  seems  Co  be  the  first  attempt  to 
investigate  Che  evolution  of  a  tropical  cloud  cluster  with  the  inclusion  of  the 
A-S  scheme  in  the  primitive  prognostic  equations.  We  are  pleased  not  only  that 
the  simulated  results  bear  considerable  resemblance  to  the  observations  but 
also  chat  this  study  reveals  some  intrinsic  nature  of  the  A-S  scheme  as  the 
convective  system  evolves. 

As  far  as  future  research  is  concerned,  we  divide  our  intention  into  four 
steps.  The  first  one  is  to  apply  the  same  semi-prognostic  method  to  a 
mldlaCiCude  data  set  to  examine  whether  or  not,  or  to  what  extent,  the  quasi- 
equilibrium  hypothesis  works  in  the  extratropical  atmosphere.  The  second  step, 
which  certainly  depends  upon  the  performance  of  the  first  one,  is  to  use  Che 
same  mesoscale  model  Co  simulate  extratropical  mesoscale  convective  systems. 


,*at  ,1 


-194- 

The  third  step  Is  to  apply  the  framework  of  the  mesoscale  model  but  with 
coarser  grids  to  easterly  wave  disturbances  to  investigate  the  ability  of  the 
model  in  reproducing  the  tropical  large-scale  phenomena.  The  fourth  step  is  to 
incorporate  the  A-S  scheme  into  a  regional  prediction  model  or  general 
circulation  model  to  test  the  ability  of  the  scheme  in  simulating  the  observed 
synoptic-scale  or  global  distributions  of  cloud  effects  and  precipitation.  As 
should  be  expected,  the  accomplishment  of  this  task  is  also  contingent  upon 
other  physical  processes,  especially  the  boundary  layer  parameterization,  which 
must  be  able  to  resolve  well  the  coupling  processes  between  the  cloud  layer  and 
boundary  layer  (Suarez  et  al. ,  1983). 


APPENDIX  A 


CLOUD  EQUATIONS 


A.l  The  cloud  model  for  cumulus  sub-ensembles; 


The  one-dlmenslonal  steady-state  entraining  cloud  model  used  In  the  A-S 


scheme  can  be  represented  by  the  following  budget  equations: 


3n(z,X) 


static  energy: 


Xn(2,x) 


(A.l) 


[n(z,x)s-(z,x)l 


moist  static  energy: 


^  h(z,X)h.(z,X)] 
3z  c 


s(z)  +  Ln(z,X)c(z,X) 
3z 


3ti(z,X) 


(A.2) 


(A. 3) 


water  content: 


—  ((q-(z,X)  +  £(z,X))] 


az 


5(z)  -  n(z,X)r(z,X) 

o  Z 


(A.4) 


where  c(z,X)  Is  the  condensation  rate,  serving  as  a  source  for  the  static 
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energy  and  a  sink  for  water  vapor  and  r(2,X)  is  the  rainwater  converted  from 
the  liquid  water  i(z,X)  by: 


r(z,X)  =  CgiCz.X) 


(A. 5) 


where  Cq  is  the  conversion  coefficient  and  is  specified  by  the  value  of 
2  X  10”^  cm“^.  The  total  rainwater  per  unit  height  (R)  can  be  calculated  by: 


,Xn(z) 

R(z)  =  /  “  m„(X)n(z,X)r(z,X)dX 

J  0  “ 


Solving  (A.l),  we  obtain  the  normalized  mass  flux: 


(A. 6) 


n(z,X)  =  exp{X(z  -  Zd)} 


With  the  use  of  (A.l),  we  can  solve  h^.  from  (A.3)i 


hj.(z,X)  =  y  (hj,(zg,X)  +1/2  n(z' ,X)H(z')dz' ] 

n  V  z  j  A  j  B 


(A. 7) 


(A. 8) 


Assuming  that  the  air  Is  saturated  in  the  clouds  and  the  pressure  difference 
between  the  cloud  and  environment  is  neglected,  we  obtain  two  approximations 
(Arakawa  and  Schubert,  1974); 


S(,(z,X)  -  s(z)  =  (h^(z,X)  -  h  ) 


(A. 9) 


qj.(z,X)  -  q  (z) 


(A. 10) 


f  (h-(z,X)  -  h  ) 

1+r  L  ^ 


C  '  -^P 

P  3T 


(A. 11) 


substituting  (A. 8)  Into  (A. 9)  and  (A.  10),  we  can  respectively  obtain  s^.  and  q^,, 
and  subsequently  substituting  q^  Into  (A. 4),  Jt(z,X)  can  be  solved.  With  the 
known  a(z,X),  we  can  determine  J,(z)  in  (2.7)  and  (2.8)  by: 

l(z)  »  i(z,Xjj(z))  (A. 12) 

where  Xq(z)  represents  the  fractional  entrainment  rate  of  the  sub-ensemble 
which  has  the  detrainment  level  at  z.  The  schematic  diagram  of  the  vertical 
structure  of  a  sub-ensemble  and  the  arrangement  of  the  cloud  variables  is  shown 
in  Fig.  A.l. 

A. 2  The  determination  of  X: 

Figure  A. 2  shows  the  schematic  diagram  of  the  distribution  of  the  cloud 
spectrum  for  a  cumulus  ensemble.  To  determine  X  for  each  sub-ensemble,  we 
first  prescribe  the  cloud  top  nonbuoyancy  condition  by: 

hj.(l)  -  h*(l)  (A. 13) 
where  h  (1)  represents  the  h^  at  the  top  of  the  1th  sub-ensemble  which  detrains 


Fig.  A.l:  The  schematic  diagram  of  the  vertical  structure  of  the  sub-ensemble 
cloud  i  where  n,  q^,  and  are  located  at  half-integer  levels  and  condensation 
rate  (c)  and  vertical  convergence  of  heat  and^moisture  fjuxes  (F)  located  at 
integer  levels.  The  cloud  top  conditions  D(i),  h^(l),  and  1(1)  are  located  at 
level  i.  The  thin  arrows  denote  the  entrainment  from  the  environment.  The 
vertical  thick  arrow  denotes  the  cloud  base  mass  flux  originating  from  the  PBL. 
The  horizontal  thick  arrow  denotes  the  detrainment  of  cloud  air. 


k 


Fig.  A.2;  The  schematic  diagram  of  the  distribution  of  cloud  spectrum  for  a 
cumulus  ensemble  In  which  the  entrainment  rates  have  the  relation  as 

^16  ^  ^15^ . also  Indicated  that  the  smaller  X  has  the 

larger  cross  section  of  the  cloud. 


at  the  level  1  and  h  (1)  is  the  large-scale  saturated  moist  static  energy  at 
the  level  1.  Then,  using  (A. 13)  and  (A. 8)  and  inserting  the  cloud  base 
condition 


b-(zD,X)  *  hv 


(A. 14) 


Into  (A. 8),  we  can  iteratively  determine  the  fractional  entrainment  rate  for 
the  ith  sub-ensemble  (X^)  with  the  use  of  (A. 7).  Thus,  it  is  evident  that  the 
value  of  h^  in  (A. 14)  is  very  crucial  in  determining  the  properties  of 
sub-ensembles . 


A. 3  The  cloud  model  for  cumulus  downdrafts: 


Since  we  model  the  cumulus  downdraft  practically  as  an  "inverted"  updraft 
(see  Fig.  A. 3),  the  cloud  equations  for  downdrafts  can  be  easily  deduced  from 
those  for  updrafts  presented  in  Section  A.l: 


static  energy: 


“Xn  z  ,x ) 


(A. 15) 


[hjjC  z  ,x  z  ,x )  ] 


moist  static  energy: 


3Tid(z,X) 

- - -  s(z)  -  nd(2,X)E(z,X) 

0  z 


(A. 16) 


Fig.  A. 3:  The  schemaclc  diagram  of  the  vertical  structure  of  a  downdraft 
element.  The  details  of  the  accompanying  updraft  have  been  Illustrated  In 
Fig.  A.l.  e„  denotes  the  evaporation  rate  In  the  downdraft. 


(A. 17) 


a  3i1  \  2  fX  ) 

^  [nd(z.X)h^d(*,X)]  -  - H(z) 

water  vapor: 

^  [nd(z.x)qcd(®»^^l  ■  — - 5(2)  +  nd(2.X)E(z,x)  (A. 18) 

where  E(z,X)  is  the  evaporation  rate  from  the  falling  precipitation  and  Is  a 
sink  for  static  energy  and  a  source  for  water  vapor.  It  Is  seen  that  a 
negative  sign  Is  In  front  of  X  In  (A. 15),  which  Illustrates  the  nature  of  the 
"Inverted  plume."  It  Is  also  seen  that  there  Is  no  liquid  water  term  (1)  In 
(A. 18)  because  no  condensation  takes  place  In  the  downdrafts.  Since  the 
downdrafts  are  also  assumed  to  be  saturated,  (A. 9)  and  (A. 10)  are  still  valid 
In  determining  s  ^  and  from  hj.j. 


APPENDIX  B 


DERIVATIONS  OF  K(X,X')  and  F(X) 


Starting  from  the  definition  of  cloud  work  function: 


A(X)  -  / 


Zd(X) 


— - - n(z,X)[8  (z,X)  -  i(z)]dz 

CpT(z) 


and  taking  the  time  derivative  of  (B.l),  we  have: 


3A(X) 

at 


J  ^  — - - n(z.X)  ^  [s  (z,x)  -  i(z)]dz 

CpT(z)  3t  c 


Using  (A. 9),  we  obtain: 


aA(X) 

at 


J^B 


(hj.(z,X)  -  h(z)ldz 


Using  (A. 8)  and  (A. 14)  yields: 


aA(X) 

at 


Zd(X)  g  1  3  fh  +  1 

^“'B  CpT(z)  1+y(z)  3t  M 


/'  n(z',X)h(z)dz' 

'b 


-  n(z.X)h*(z)]dz 


By  the  relation  of 


where  y  has  been  defined  in  (A. 11),  (B.4)  can  be  expressed  as: 


at  ■’ CpT(z)  i-hrC*) 


at 


at 


-  n(2,X)(l  +  y(z))  ^^Idz 


at 


Reversing  the  order  of  integration  of  the  double  integral  in  (B.6),  we  have 


at  a t  J  Zb  p 


dz 


ah(z) 


at 


dz 


-/ 


Zd(A )  g 


^B  CpT(z) 


n(z,x) 


as(z) 


at 


dz 


Now,  by  defining: 


c„T(z)  l-n(z) 


=  p(z)6(z) 


(B.12) 


Substituting  Mj.(z)  and  D(z)  from  (2.10)  and  (2.16)  into  (B.12),  we  obtain 

-  Fj^(X)  +  F^(X)  +  /q“**[K^(X,A')  +  KD(X,X')]mB(X')dX' 

3  t 

where 

Fc(X)  “  /z°^  ^  n(z,X){xb(z,X)[-w  -  t*7h  +  Qj^] 


-  (i(z)(-  w  —  -  ^*78  +  Qj^)}  p(z)dz 
9  z 


9  h»i  Zn(X  ) 


K^(X,X')  -  n(z.X)n(z.X'){Xb(z,X)  ^  -  a(z)  —}  dz 


U-momentum  equation: 


„n+l  „n-l 

“1+1/2, k  ■  “1+1/2, k 

2At 


^  ^“1+1/2, k  “1+3/2, k^  2 


,.“1+1/2, k  “1-1/2, k^2 

'■ - 2 -  /^*l+l/2 


,  r,''l,k+l/2  ''l+l,k+l/2s  ,“i+l/2,k  “i+l/2,k+l 

l'’k+l/2  U - 5 - )  I - 2 - 


r(."l,k-l/2  ''i+l,k-l/2w “1+1/2, k  +  “1+1/2, k-K, 

“  Pk-1/2  U - 2 - - 2 - 


^''1+1/2, k  ■  ‘^p®k^’'i+l,k  ■  ’'l,k^''^*l+l/2 


v-monentum  equation: 


where  n  denotes  the  time  step,  i  and  k  respectively  denote  the 
vertical  grid  points,  and 

^*i+l/2  “  ^^*i 

^^k+1/2  “  ^^^k 


horizontal  and 


(C.5) 


(C.6) 


Note  that  the  time-independent  imposed  uq  and  Wq  are  not  included  in  the  above 
equations  for  simplicity  and  the  moisture  equation  is  also  not  included  since 
it  has  the  same  form  as  (C.A). 
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